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ABSTRACT 


The  analytical  model  for  the  erosion  of  a  brittle  mate- 
rial  by  multiple  particle  impacts^  is  revised  in  accordance 
with  microscopic  observations  of  the  erosion  damage  due  to  the 
direct  impact  of  glass  beads.  Procedures  are  described  for 
obtaining  the  quantitative  data  required  to  make  the  model  a 
predictive  tool.  Erosion  tests  were  conducted  in  the  AFML- 
Fell  erosion  facility  using  70  and  290  micron  diameter  glass 
beads  and  uniformly-graded  quartz  sand  impacting  flat-faced 
specimens  of  borosilicate  glass  directly  and  at  an  oblique 
angle.  Excellent  agreement  was  found  between  the  average  ring, 
fracture  diameter  and  the  diameter  of  the  contact  area  pre¬ 
dicted  by  the  Hertzian  theory  of  imract  for  the  glass  bead 
impacts.  The  damage  mechanisms  for  glass  plates  Impacted  by 
sand  particles  and  liquid  drops  were  investigated  through 
microscopic  examination  of  the  progressively-eroded  surfaces. 

Weight-loss  data  are  available  for  ultra-high  molecular 
weight  polyethylene*  nickel-coated  glass-fiber  reinforced 
epoxy*  and  graphite-fiber  reinforced  epoxy  specimens  exposed 
to  rain  and  sand  erosion  environments.  The  polyethylene 
exposed  to  a  rain  environment  eroded  in  a  preferential  manner 
as  hypothesized  in  the  analytical  model  being  developed  for 
glass  plates  based  on  the  concept  of  pit  nucleation  and  growth. 
Cutting  wear  was  the  primary  erosion  mechanism  for  sand  ero¬ 
sion  of  the  UHMW  polyethylene.  Nickel-corted  composite  speci¬ 
mens  with  1  and  6  mil  coatings  exposed  to  the  sand  erosion 
environment  and  the  3  mil  coating  in  the  rain  environment 
exhibited  buckling  of  the  coating.  Microscopic  examination 
of  these  specimens  did  not  provide  Insight  into  the  critical 
conditions  required  to  produce  buckling  of  the  coatings*  The 
erosion  mechanisms  associated  with  the  intermediate  to  long 
time  rain  and  sand  erosion  of  the  graphite-fiber  reinforced 
epoxy  specimens  have  been  Identified,  A  micromechanical  model 


for  material  removal  due  to  liquid  droplet  impacts  is  proposed 
using  simple  beam  theory.  The  critical  material  parameters 
In  this  case  are  values  of  the  mode  I  fracture  toughness  evalu¬ 
ated  for  an  initial  crack  lying  parallel  and  transverse  to 
the  graphite  fibers.  Examination  of  the  graphite- fiber  rein¬ 
forced  epoxy  specimens  exposed  to  a  sand  environment  reveals 
that  both  cutting  and  direct  deformation  erosion  are  equally 
effective  in  removing  fragments  of  material  from  the  surface 
which  at  most  contain  a  few  short  segments  of  fiber  and  the 
surrounding  matrix. 
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1.0  INTRODUCTION 


The  initiation  stage  of  the  general  erosion  process  for 
brittle  and  semibrittle  materials  can  be  idealized  and  des¬ 
cribed  in  qualitative  terms  as  follows.  Particles  impact  a 
solid  material  surface.  The  bulk  material  has  a  certain 
microstructure  and  a  certain  flaw  distribution.  The  form  of 
the  microstructure  may  be  the  grains  in  a  polycrystalline 
solid,  crystalline  phases  in  a  semicrystalline  polymer,  or  a 
disperse  phase  in  a  composite  such  as  the  fibers  in  a  fiber- 
reinforced  composite.  The  flaw  distribution  may  be  pre^ipi - 
tates  in  ductile  metals,  microscopic  surface  cracks  in 
inorganic  glasses,  or  delaminated  regions  at  the  fiber/matrix 
interface  in  composites.  The  localized  loadings  are  not  sus¬ 
tained  but  generate  stress  pulses  which  propagate  through  the 
material.  These  stress  pulses  can  induce  minute  flaws  in  the 
material  if  they  are  not  already  present  through  interactions 
with  the  microstructural  features  of  the  material,  or  they 
can  extend  existing  flaws  to  a  point  where  they  become  unstable. 
This  process  continues  as  a  random  function  of  time  until 
localized  flaw  distributions  reach  a  critical  stage  whereby 
pieces  of  material  are  freed  from  the  surface.  The  efficiency 
with  which  material  removal  is  initiated  will  depend  on  the 
magnitude,  duration,  and  frequency  of  the  localised  loadings 
and  the  material's  response  to  these  loadings. 

In  order  to  understand  the  erosion  process  in  these  terms 
and  to  obtain  analytical  evaluations  for  the  <  rosion  of  speo-i* 
fit  materials,  three  major  areas  of  investigat  ion  must  be  con-., 
sidered: 

•Characterization  of  the  localized,  time-dependent 
pressure  distribution  applied  to  a  material  surface  ..  . 
by  an  impacting  particle,. 


•Characterization  of  the  transient  stress  distribution 
within  the  body  due  to  the  impacting  particles;  pro¬ 
pagation  of  stress  and  shock  waves. 

•Characterization  cf  the  dynamic  fracture  behavior  of 
materials  based  on  the  interaction  of  propagating 
stress  waves  with  the  material's  microstructure  and 
distributed  defects. 

The  complex  conditions  described  here  are  only  an  indication 
of  the  overall  problem.  The  objective  in  our  modeling  studies 
has  been  to  find  justification  for  reasonable  simplifications 
in  the  genrral  analyses  required.  Evidence  has  already  been 
provided  for  inorganic  glass  plates  impacted  by  spherical 
glass  beads  that  certain  results  from  the  Hertzian  (quasi- 
static)  theory  of  impact  are  applicable  to  a  problem  which 
involves  highly  transient  stress  distributions.  In  the  area 
of  fracture  mechanics  we  are  attempting  to  delimit  those  con¬ 
ditions  under  which  the  statically  determined  fracture,  data . 
can  be  applied  to  dynamic  fracture,  problems;.  Results  of  this 
type  provide  a  rational  basis  for  expecting  correlations 
between  erosion  behavior  and  the  mechanical  properties  of. 
particular  materials  which  are  determined  from  micromechanical, 
analyses  based  on  microscopic  examination  of  the  erosion  damage 
at  various ■ stages  of. the  overall  erosion  process,  failure  to 
examine  the  erosion  process  at  the  level,  we  are  suggesting 
lias  resulted  in  the  general .  sparcit y  of  correlations  between  • 
erosion  and  mechanical  behavior  as  exemplified  by  the  negat i v- 
results  of  numerous  investigat ora..  ;  While  the  research  report  ed 
here  has;  not  completely  solved  this  problem,  considerable  htOVi* 
way  has  been  made  in  understanding  and  analyzing  the  erosion 
behavior  of  moderately  complex  material  systems. 


2.0  SUMMARY 


A  portion  of  the  research  reported  here  is  a  continua¬ 
tion  of  research  described  by  Adler  and  Sha.  ^  The  micro¬ 
scopic  observations  of  70  and  290  micron  glass  beads  impacting 
glass  plates  have  provided  a  complete  picture  of  the  erosion 
process  for  the  entire  range  of  particle  sizes  which  are  of 
interest  in  dust,  sand,  and  ice  particle  erosion  at  subsonic 
velocities. 

The  general  formulation  of  an  analytical  model  describ¬ 
ing  the  erosion  of  brittle  materials^1  \  has  been  made  more 
explicit  in  accordance  with  the  experimental  findings  per¬ 
taining  to  the  erosion  of  glass  plates.  The  present  research 
establishes  experimental  and  analytical  procedures  for  the 
attainment  of  quantitative  data  to  specify  the  functions 
required  in  the  revised  equations  governing  the  erosion  pro¬ 
cess.  The  results  of  .the  erosion  tests  using  boros! licate 
glass  plates  (horning  Glass  Works  designation  Pyrex  7740)  for 
the  erosion  targets  indicated  that  the  diameters  of  the  >-ing 
fractures  produced  on  the  specimen  surfaces  by  both  the  70 
and  290  micron  glass  heads  impact ing  over  a  velocity  range 
cf  from  200  to  kl\0  fp«r  were  in  excellent  agreement  with  the 
diameters  of  the  maximum  contact  areas  calculated  from  the 
Hertzian  theory  of  impact.  This  fact  implies  that  the  tran¬ 
sient  nature  of  the  bead  impacts  is  of  secondary  importance 
in  the  specification  of  the  Information  required  to  describe 
the  general  erosion  process. 

The  weight-loss  data  obtained  for  oblique  impacts  by 
spherical  glass  beads  at  the  same  normal  velocity  used  in  the 
case  of  direct  impacts  on  glass  plates  indicate  that  under 
nearly  equivalent  conditions  the  magnitude  of  the  weight  loss 
at  the  same  exposure  time  for  oblique  impacts  was  two  to  three 
times  greater  than  the  corresponding  value  for  direct  impacts. 
The  ring  fracture  diameters  were  found  to  be  approximately 
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equal  for  both  impact  conditions,  however  the  penetration 
depths  of  the  cone  fractures  were  noticeably  less  for  the 
oblique  impacts  especially  as  the  impact  velocity  increased. 

It  is  concluded,  for  a  reason  yet  to  be  discovered,  that  the 
probability  for  producing  ring ' fractures  for  the  oblique 
impacts  was  considerably  higher  than  for  direct  impacts. 

Exploratory  investigations  for  the  direct  and  oblique 
impact  of  quartz  sand  particles  on  glass  plates  were  initiated 
For  both  impact  conditions  the  normal  component  of  the  par¬ 
ticle  impact  velocity  was  held  constant.  It  was  found  that 
the  steady-state  erosion  rates  were  nearly  identical.  The 
slightly  greater  values  obtained  for  oblique  impacts  may  be 
due  to  the  fact  that  the  specimen  receives  approximately  30 
percent  more  particle  impacts  per  revolution  in  this  case. 
Microscopic  examination  of  the  damage  occurring  during  the 
early  a  - ages  of;  the  erosion  process  did  not  reveal  any  signi¬ 
ficant:  differences  between  the  direct  and  oblique  impacts. 

As  would  be  expected  the  erosion  rates  for  the  If regularly- 
shaped  particles  are  considerably  higher  than  those  for  glass 
beads  of  an  equivalent  dimension.  The  basic  requirements  for 
characterizing  the  particle  shapes  in  a  distribution  of 
irregularly-shaped  particles  are  described  in  conjunction  with 
the  extension  of  the  analytical  mod^l  to  erosion  of  brittle 
materials  by  sand  particles, 

Itie  initial  damage  observed  on  the  flat  surface  of  a 
glass  specimen  arising  from  liquid  droplet  impacts  is  due 
primarily  to  the  transient  pressure  d 1st ributlon  locally 
applied  to  >he  surface.  Plane  glass  surfaces  offer  the  advan¬ 
tage  of  an  inherently  simple  mechanical  response  to  i he . 
appl >d  loadings  which  closely  approximates  that  considered 
In  most  of  the  mat heretical  analyses  of  stress  wave  propaga¬ 
tion  in  an  elastic  half-space.  A  procedure  Is  outlined  for 
estimating  the  magritude  of  the  time-  and  spatially-depend^nt 
pressure  pulse  applied  to  a  deformable  r***r*co  during  the 
Initial  stages  of  a  spherical  liquid  drop  striking  a  material 
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surface  at  low  to  moderate  velocities.  The  interaction 
between  the  deforming  solid  surface  and  the  liquid  drop  is 
taken  into  account  and  the  transient  stress  distribution  with¬ 
in  the  solid  target  can  be  evaluated.  The  erosion  tests  on 
•glass  plates  provided  insights  into  the  fracture  patterns 
which  developed  due  to  droplet  impacts  at  300  to  500  fps. 

Direct  correlations  with  the  numerical  evaluation  of  the  trans¬ 
ient  stress  distributions  which  could  be  made  on  the  basis  of 
the  above  analysis  are  not  available,  since  the  acquisition 
of  this  information  requires  fairly  extensive  computer  pro¬ 
gramming  which  will  be  initiated  in. the  future.  The  erosion 
tests  established  that  material  removal  in  a  rain  environment 
was  a  preferential  process  of  pit  nuoleat ion  and  growth,  ; 
however  weight- loss  data  are  not  available  due  to  cracking 
of, the  specimens  which  occurred  before  a  sufficient  number 
of  data  points  could  be  obtained.  This  experimental  obser¬ 
vation  and  the  transient  stress  wavs  analysis  described  above 
provide  a  basis  for  extending  the  analytical  model  developed 
f  ,r  solid  particle  impacts  to  liquid  drop  impacts. 

Height- loss  data  are  available  for  ultra-high  molecular 
weight  (tJHMWl  polyethylene,  nickel-coated  glass-fiber  epoxy* 
and  graphite- fiber  reinforced  epoxy  specimens  exposed  to  rain 
and  sand  erosion .Environments .  A  number  of  specimens  eroded 
to  various  stages  of  the  overall  erosion  process  are  also 
available  for  direct  microscopic  identification  of  the  opera¬ 
tive  erosion  mechanisms  supplementing  the  microscopic  inves¬ 
tigations  which  were  conducted  while  the  specimens  were  being 
Incrementally  eroded  in  the  AFML- Beil  ratal ing-arm  facility. 

The  rain  erosion  characteristics  of  IMMW  polyethylene 
exhibited  a  strong  tendency  toward  a  process  of  pit  nuclea? ion 
and  growth.  Over  a  velocity  range  of  from  500  to  112C  fps 
the  time  scale  for  which  a  particular  erosion  mechanism  is 
operative  in  these  specimens  is  only  contracted  or  expanded 
as  the  impact  velocity  is  increased  or  decreased  but  no  -hange 


is  evident  in  the  sequence  of  operative  erosion  mechanisms. 
Material  removal  from  aerodynamically- shaped  specimens  by 
sand  particles  is  due  to  the  cutting  action  of  the  particles: 
more  material  is  removed  at  locations  on  th  specimen’s  sur~ 
face  away  from  the  leading  edge  than  at  the  leading  edge. 

A  steady-state  erosion  rate  is  quickly  established  in  this 
case. 

The  nickel-coated  composites  exposed  to  a  rain  environ¬ 
ment  developed  a  dense  array  of  cracks  in  the  coating  which 
rarely  penetrated  the  full  thickness  of  the  coating.  The 
influence  of  the  coating  thickness  on  erosion  behuvior  was 
investigated.  A  specimen  with  a  3  mil  coating  failed  by 
local iced  buckling  of  the  coating  followed  by  catastrophic 
failure  of  the  entire  buckled  area.  At  this  time  a  few  ero¬ 
sion  pits  could  be  found  oh  the  surface  which  almost  penetra¬ 
ted  the  coating  thickness.  A  specimen  with  a  6  mil  coating, 
did  not  exhibit  the  buckling  instability,  but  rather  the 
Coating  was  completely  eroded  in  a  localized  region  on  the 
leading  edge  of  the  specimen  at  an  exposure  time  which  was 
seven  times  longer  than  that  at  which  buckling  occurred  for 
the  3  mil  coating*  No  measurable  weight  loss  was  observed  in 
either  ease.  Specimens  with  both  the  3  and  6  mil  nickel  coat¬ 
ings  were  exposed  to  a  sand  environment  at  ?30  fps.  Both 
coating-  thicknesses  failed  by  buckling  of  the  coat ings  which 
extended  over  approximately  one-half  the  length  of  the  speci¬ 
men.  Measurable  weight  losses  were  recorded  in  both  cases. 

The  -rate  of  material  removal  was  consent  after  a  short incu¬ 
bation  period  and  essentially  identical-  for  both  specimens. 

The  6  mil  coating  lasted  three  times  longer  than  -he  3  mil 
coating  in  the  sand  environment,  and  the  specimen  with  the 
6  mil  coating  had  a  weight,  loss  approximately  three  times  . 
the  weight  !o3$  for  the  3  mil  coating  at  the  tine  of  failure. 
Microscopic  examination  of  the  nickel-coated  compos  it «?  speci¬ 
mens  did  not  provide  Insight  Into  the  critical  conditions 
required  to  produce  buckling  of  the  coating. 


The  erosion  mechanisms  associated  with  the  intermediate 
to  long  time  rain  and  sand  erosion  of  the  graphite -fiber 
reinforced  epoxy  specimens  have  been  identified.  For  a  rain 
environment  a  micromechanical  model  is  presented  for  material 
removal  based  on  simple  beam  theory.  According  to  this  model 
the  important  material  characteristics  involved  in  the  erosion 
behavior  of  the  composite  are  its  mode  I  fracture  toughness 
values  evaluated  under  dynamic  conditions  for  an  initial  crack 
parallel  to  the  fiber  reinforcement  and  transverse  to  the 
fiber  reinforcement.  The  sand  erosion  of  the  graphite-fiber 
reinforced  epoxy  specimens  occurs  through  fine- scale  cracking 
of  both  the  fiber  and  matrix.  The  eroded  fragments  would  con¬ 
tain  at  most  a  few  short  segments  of  *  he  fibers  and  the  sur¬ 
rounding  matrix.  A  steady-state  erosion  rate  is  quickly 
established  at  an  impact  velocity  of  <30  fps  after  a  very  short 
incubation  period.  On  the  basis  of  the  uniformity  of  the 
eroded  surface  of  the  a<  rcdynamioaUyr-shaped  specimens,  direct 
deformation  and  wear  deformation  ar*  equally  effective  in 
removing  material  from  the  surface  of  this  composite. 


3.0  EROSION  OF  GLASS  SPECIMENS 

Erosion  data  were  obtained  for  borosilicate  glass  speci¬ 
mens  (Corning  Glass  Works  designation  Pyrex  7740)  subjected 
to  the  direct  and  oblique  impact  of  70-micron  glass  beads, 
290-micron  glass  beads,  and  quartz  sand  with  a  mean  diameter 
of  approximately  90  microns.  The  erosion  tests  conducted  are 
recorded  in  Table  I,  Due  to  the  scope  of  the  present  program 
it  was  only  possible  to  run  one  or- at  most  two  specimens  under 
the  same  test  conditions,  so  the  findings  reported  here  are 
only  indications  of  the  general  features  of  the  erosion  process. 
These  results  are  adequate  however  for  providing  guidance  in 
extending  the  solid-particle  modeling  studies  to  more  general 
erosion  conditions. 

3*1  General  Considerations 

Glass  specimens  were  used  in  the  erosion  tests.  The 
mechanical  properties  of  Pyrex,  fused  silica,  and  soda  lime 
glass  are  listed  In  Table  II.  Table  III  indicates  the  magni¬ 
tude  of  the  relevant  indentation  parameters  based  on  the 
Hertzian  theory  of  impact.  To  slide  rule  accuracy  the  hange 
in  the  magnitude  of  the  impact  parameters  is  negligible  for 
Pyrex  impacted  by  soda  lime  g lass  beads  compared  with  fused 
silica  impacted  by  soda  lime  glass  beads*  Microscopic  inves¬ 
tigation  of  the  general  features  of  the  erosion  process  in 
each  case  show  them  to  be  identical. 

The  glass  specimens  were  rectangular  parallelepipeds: 
3.25  ih«  x  0,875  in.  x  0.5  in.  In  order  to  obtain  a  uniform 
distribution  of  particle  impacts  on  the  face  of  the  specimen 
and  to  eliminate  edge  effects,  thin  metal  plates  were  inserted 
in  the  specimen  holder  so  that  only  the  central  portion  of  the 
specimen's  surface  is  exposed  to  the  particles  dropping  from 
the  nosale  at  the  end  of  the  blast  tube  in  the  AFML-Sell  ero¬ 
sion  facility.  The  specimen  configurations  for  direct  impacts 
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TABLE  II 


MECHANICAL  PROPERTIES 
OP 

GLASS  BEADS  AND  SPECIMENS 


Soda  Lime 
Glass 

Fused  Silica 

Pyrex 

Density 

(lb. /in.3) 

0.0866 

0:.  0795 

0.0805 

Poisson1 

1  s  Ratio 

0.22 

0.17 

0.22- 

Young's 

Modulus  (lb. /in.2) 

10.0  x  106 

10.68  x  106 

9.8  x  lO1 

TABLE  III 

HERTZIAN  IMPACT  PARAMETERS 


Equivalent  Static 
Load,  P(lb. ) 

Radius  of  Contact 
Circle,  a (in.) 

Duration  of 
Contact, T (sec) 

70  Micron 

Glass 

Beads 

200  fps 

0.215 

0.350 

X 

10“  3 

1.19  x 

10"  7 

300  fps 

0.355 

0.415 

X 

10“  3 

1.10  x 

10"  7 

400  fps 

0.504 

0.465 

X 

10“  3 

1.04  x 

10"7 

500  fps 

0.653  . 

0.506 

X 

10"  3 

0.99  X 

10"  7 

290  Micron. 

Glass 

'  • 

Beads 

200  fps 

3.56 

1.430 

X 

10"  3 

4.76  X 

10"  7 

300  fps 

5.88 

1.692 

X 

10“  3 

4.40  x 

10"  7 

'400  fps 

8.35 

0 

0 

CR 

H 

X 

10"3 

4.15  x 

10”7 

500  fps 

10.80 

2.078 

X 

10"  3 

3.97  x 

O 

1 

-<] 
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and- oblique  impacts  when  the  normal  to  the  specimen's  surface 
is  inclined  at  an  angle  of  50°  to  the  plane  of  the  rotating 
arm  are  indicated  in  Fig.  1.  For  direct  impacts  the  exposed 
area  on  the  face  of  the  specimen  is  0.82  in.  (1.05  in.  x 
0.78  in. )  and  for  the  oblique  impacts  it  is  1.8l  in . ^  (0.59  in. 
x  3,0?  in.). 

A  sieve  analysis  was  made  of  the  70-micron  glass 
beads  and  a  quartz  sand  supplied  by  the  American  Graded  Sand 
Company,  Paterson,  New  Jersey,  which  was  to  be  narrowly  dis¬ 
tributed  around,  a  mean  particle  size  of  70  microns.  The  sieve 
sizes  used  are  listed  in  Table  IV,  and  a  comparison  of  the 
size  distribution  for  the  sand  and  glass  beads  is  shown  in 
Fig.  2.  It  is  seen  that  the  diameters  of  the  glass  beads  are 
quite  narrowly  distributed  around  e.  mean  diameter  of  70  microns. 
The  mean  maximum  lateral  dimension  of  the  quartz  sand  particles 
is  estimated  to  be  92  microns;  the  size  distribution  is  not 
very  narrow,  but  fairly  smoothly  spread  over  a  100-micron  range 
above  63  microns.  The  angularity  of  the  sand  particles  can 
be  seen  in  Fig.  3*  As  a  matter  of  convenience  in  writing,  the 
quartz  sand  will  be  referred  to  as  the  90-micron  sand. 

Modification  of  the  formula  derived  earlier^1)  for 
estimating  the  number  of  particles  impacting  the  surface  of 
the  specimen  is  required  to  Include  the  reduction  in  the  sur¬ 
face  area  of  the  specimen  exposed  to  particle  Impacts  and  the 
inclusion  of  particles  striking  the  surface  of  the  specimen 
obliquely.  The  final  result  is  simply 


N 

720s/2gK 


d'cosfl 


where  is  the  number  of  particle  impacts  the  specimen 
experiences  per  revolution  of  the  rotating  arm' 

N  is  the  number  of  particles  released  per  minute 
sfSg H  3  21.233  ft.  for  the  AFML-Bell  erosion  facility 
1  is  the  lateral  dimension  of  the  nozzle  in  inches 
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IGURS  1  liozzle  dimensions  and  experimental  arrangement  for  direct 

and  oblique  solid  particle  impacts. 


TABLE  IV 


SIEVE  SIZES  USED  IN  THE  DETERMINATION  OF 
PARTICLE  SIZE  DISTRIBUTION  FOR  GLASS  BEADS 
AND  QUARTZ  SAND 


Sieve  Number 

Tyler  Equivalent  Mesh 

Mesh  Opening 

Microns 

Inches 

80 

177 

0.0070 

.140 

-  .105  .  .... 

.0041 

170 

90 

.0035 

230 

63 

.0025 

270 

-53 

.0021 

.  325 

44 

.0017 

90-aicrnn  ■sand 


I'  is  the  lateral  dimension  of  the  exposed  area  on 
the  specimen's  surface  in  inches 
d'  is  the  transverse  dimension  of  the  exposed  area  on 
the  specimen's  surface  in  inches 
0  is  the  angle  of  inclination  between  the  normal  to 
the  specimen's  surface  and  the  plane  of  the  rotating 
arm. 

When  0  =  90°,  -  0  according  to  Eq.  (1).  This  result  is 

consistent  with  the  implicit  assumption  that  the  vertical 
velocity  of  the  particles  is  negligible  with  respect  to  the 
velocity  of  the  specimen. 

For  the  experimental  arrangement  illustrated  in 
Fig.  1,  when©  =  0°,  £  =  1.375  in.,  £'  =  1.05  in.,  and  d'  = 
0.78  in.  and  when  0=  50°,  £  -  1.50  in.,  £'  =  0.59  in.,  and  > 
d'  =  3*07  in.  It  has  been  assumed  that  the  lateral  dimension 
of  the  nozzle  is  representative  of  the  equivalent  uniform 
particle  distribution  across  the  specimen  for  the  entire  vol¬ 
ume  of  particles  dropping  through  the  path  of  the  specimen. 
Admittedly  this  approximation  does  not  properly  account  for 
the  observed  lateral  flaring  of  the  particle  flow  as  i+  leaves 
the  nozzle,  however  the  uniformity  of  the  flow  near  the  center 
of  the  specimen  does  not  seem  to  be  affected  by  the  decrease 
in  the  particle  flow  density  at  the  ends  of  the  specimen.  The 
results  of  Eq.  (1)  should  therefore  provide  a  reasonable  esti¬ 
mate  of  the  number  of  bead  impacts  per  revolution  for  the 
experimental  arrangement  used  in  the  erosion  tests. 

For  the  glass  beads  a  reasonable  estimate  can  be 
made  for  the  number  of  impacts  the  specimen  would  experience 
for  a  given  bead  flow  rate  and  velocity.  Further  assumptions 
are  required  to  provide  a  similar  estimate  for  the  sand  par¬ 
ticles.  The  simplest  approach  is  to  assume  that  the  sand 
particles  are  equivalent  to  90-micron  diameter  quartz  spheres 
and  compute  the  number  of  impacts  on  the  basis  of  the  formula 
in  Eq.  (1).  The  result  for  this  case  and  for  the  two  sizes 
of  glass  beads  are  listed  in  Table  V. 
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3.2  Direct  Impact  by  Glass  Beads 

The  direct  impact  of  spherical  glass  beads  on  glass 
plates  provides  idealized  conditions  for  a  combined  analytical 
and  experimental  study  of  the  erosion  process  and  its  relation 
to  the  mechanical  behavior  of  the  colliding  bodies  and  impact 
parameters.  The  mechanical  properties  of  the  target  material 
and  impacting  particles  are  those  of  a  brittle  elastic  solid. 
The  simplified  geometry  of  the  impacting  particles  allows  the 
Introduction  of  concepts  and  results  from  the  Hertzian  theory 
of  impact.  The  relations  required  from  the  Hertzian  theory 
are  discussed  by  Adler  and  Sha^^.  Additional  experimental 
and  analytical  results  which  are  essential  to  the  development 
of  the  analytical  model  proposed  by  Adler  and  Sha^^  will  be 
described. 


3.2.1  Description  of  Erosion  Process 

The  erosion  mechanisms  operative  on  flat 
surfaces  of  Pyrex  and  fused  silica  glass  specimens  impacted 
by  70  and  290  micron  glass  beads  have  been  established  through 
detailed  microscopic  investigation  of  a  number  of  progres¬ 
sively  eroded  surfaces.  The  following  patterns  emerge. 

The  basic  element  for  the  initiation  of  ero¬ 
sion  damage  in  glass  specimens  is  the  ring  fracture The 
diameter  of  the  ring  and  the  geometry  of  the  conical  fracture 
surface  extending  into  the  bulk  of  glass  plates  impacted  by 
glass  beads  Is  dependent  on  the  radius  and  impact  velocity  of 
the  glass  beads  as  well  as  the  material  properties  of  the 
colliding  bodies.  For  the  range  of  conditions  used  in  this 
program  the  conical  fracture  surfaces  are  typically  5  to  20 
microns  deep  and  on  the  order  of  1  mil  in  diameter  at  the  sur¬ 
face  of  the  glass  plate  for  the  70  micron  beads.  The  290 
micron  beads  produce  fracture  surfaces  which  are  50  to  90 
microns  deep  and  approximately  3  mils  in  diameter  at  the  sur¬ 
face  of  the  glass  plate. 
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Subsequent  Impacts  in  the  vicinity  of  the 
shallow  fracture  frustrums  produced  by  the  70  micron  beads 
initiate  chipping  of  material  from  the  surface  on  a  very  fine 
scale.  Microscopic  observations  of  the  time  evolution  of  the 
eroded  surface  indicate  that  once  a  ring  fracture  occurs  it 
is  the  nucleation  site  for  the  removal  of  material.  The 
series  of  micrographs  in  Fig.  4  and  5  show  the  change  in  the 
surface  of  the  specimen  before  a  measurable  weight  loss  is 
recorded.  The  micrographs  were  taken  in  bright  field,  so  the 
black  spots  on  the  white  field  represent  the  material  that 
has  been  removed  from  the  surface.  The  details  of  the  eroded 
surface  are  not  visible,  but  the  extent  of  the  area  removed 
is  quite  distinct.  The  eroded  areas  grow  preferentially 
around  a  ring  fracture  by  fine- scale  chipping  of  material  due 
to  intersecting  ring  fractures.  This  process  continues  to 
remove  material  from  the  surface  layer  and  expose  a  highly- 
fractured  and  irregular  subsurface  from  which  material  is 
removed  by  chipping  of  small  pieces  with  each  impact. 

In  the  case  of  the  290  micron  beads  another 
impact  in  the  vicinity  of  a  deep  conical  fracture  surface  is 
generally  not  sufficient  to  cause  material  removal  but  con¬ 
tributes  instead  to  a  more  severely  fractured  subsurface. 
Additional  impacts  in  this  region  produce  a  predictable  pat¬ 
tern  of  subsurface  fractures  at  depths  which  still  allow  a 
substantial  layer  of  overlying  material  to  remain  intact.  It 
is  only  after  a  number  of  impacts  have  occurred  in  the  vici¬ 
nity  of  the  initial  ring  fracture  that  removal  of  the  material 
exterior  to  the  developing  distribution  of  conical  frustrums 
takes  place.  The  material  contained  within  the  frustrums 
formed  by  the  conical  fracture  surfaces  is  more  resistant  to 
erosive  attack  than  the  highly  fractured  interstitial  regions, 
therefore  material  removal  from  the  interstitial  volumes  takes 
place  at  a  higher  rate  than  during  the  ring  fracture  develop¬ 
ment  phase.  The  interstitial  material  is  removed  leaving  a 
dense  distribution  of  solid  frustrums  protruding  from  the  bulk 
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(a)  1  minute  exposure 


(b)  2  minute  exposure 

Figure  4  Progressive  Erosion  of  a  Glass  Plate  by  70-Micron  Beads  at 
375  fps  11  OX  (Sheet  1  of  2) 


(a)  1  minute  exposure 


igure  fj  Erosion  of  pyrex  by  70  micron  beads 
at  375  fps  (Area  3)  110X 


specimen.  The  various  stages  in  the  materiel  ..removal  process 
are  shown  in  Pig.  6.  The  impacting  beads  chip  small  pieces 
of  material  from  these  frustrums  at  a  slower  rate  than  in  the 
interstitial  material  removal  stage.  As  the  frustrums  erode, 
more  surface  area  is  susceptible  to- bead  impacts.  Chipping 
of  the  highly- fractured  and  irregular  surface  of  the  bulk 
material  takes  place  at  a  uniform  rate  and  constitutes  the 
final  phase  of  the  erosion  process. 

The  microscopic  observations  of  the  erosion 
process  associated  with  the  two  be^d  sizes  used  in  the  ero¬ 
sion  tests  indicate  "that  the  depth  of  the  cone  fracture  sur¬ 
face  determines  whether  a  single  impact  in  the  vicinity  of  a 
pre-existing  ring  fracture  will  result  in  material  removal  or 
a  number  of  impacts  are  required  before  material  removal  is 
therefore  dependent  on  the  radius  and  velocity  of  the  impac¬ 
ting  glass  beads,  however  the  general  erosion  process  remains 
unchanged  over  the  range  of  bead  sizes  and  impact  velocities 
used  in  this  research. 

The  mechanisms  for  material  removal  described 
above  are  summarized  as  follows: 

(1)  material  removal  during  ring  fracture 
formation, 

(2)  material  removal  during  the  growth  of 
nucleated  erosion  pits, 

(3)  material  removal  due  to  chipping  of  the 
conical  frustrums  resulting  from  the 
fracture  process  in  glass  plates, 

(4)  material  removal  due  to  chipping  of  the 
bulk  material. 

Now  that  these  mechanisms  have  been  identified, 
this  information  will  be  included  in  the  mathematical  formu¬ 
lation  of  the  basic  erosion  model. 
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Figure  6 


3.2.2  Analytical  Model  for  Solid  Particle  Impacts 

The  general  erosion  process  can  be  divided 
into  three  phases.  The  total  volume  of  material  removed  is 

AV(  t  )=AV1  (t  )+AV g  ( t.  )+AV3  ( t  )  (2  ) 

where  AV-^(t)  is  the  volume  of  material  removed  at  the  time 
t  due  to  the  formation  and  growth  of  erosion 
pits  on  the  surface  of  the  specimen 
AVp(t)  is  the  volume  of  material  removed  due  to  chip¬ 
ping  of  the  conical  frustrums 
AV_(t)  is  the  volume  of  material  removed  due  to  chip¬ 
ping  of  the  bulk  material 

In  terms  of  the  previous  model,  the  general  erosion  pro¬ 
cess  can  be  described  by  dividing  the  specimen  into  two 
platelets:  a  surface  layer  defined  by  the  mean  depth  of  pene 

tration,  D,  of  the  conical  fracture  surfaces  and  the  rema4,nde 
of  the  bulk  material.  The  volume  V-^  of  the  surface  layer  Is 
now  composed  of  two  components, 

V1  =  AD  =  liHl  (AV1(t)+AV2(t))  (3) 

where  A  is  the  area  of  the  specimen  exposed  to  the  erosive 
environment . 


AV1(t)  =  r?DAA1(t)  (4) 

where  rj  is  a  shape  factor  which  will  be  described  below  and 
AA-^(t)  is  the  area  removed  at  time  t  from  the  free  surface  of 
the  first  platelet.  With  AV1(t)  defined,  the  evaluation  of 
AA1(t)  is  of  major  concern. 

The  number  of  ring  fractures  formed  on  a 
unit  area  of  the  specimen  at  any  time  t  is 


-31- 


(5) 


N(t)  =  I  (l-AA1  (t  )-N{r  )AQJn<p(  r  )dr 

where  AA^(r)  is  the  area  removed  from  a  unit  area  of  the  sur¬ 
face  layer  at  any  time  r=t  due  to  the  impacting 
particles 

Ac  is  the  area  of  a  ring  fracture 

n  is  the  number  of  impacts  /unit  area-time 

<p( t)  is  a  time-dependent  probability  that  a  ring 
fracture  will  occur. 

The  variables  are  now  expressed  in  terms  of  a  unit  of  surface 

area.  Corresponding  adjustments  should  be  made  in  Eq.  (2)  to 

(4). 


The  surface  area  of  the  specimen  exposed  to 
the  erosive  environment  changes  with  the  exposure  time,  and 
the  material  removal  due  to  the  growth  of  pits  once  they  are 
formed  must  be  subtracted  from  the  total  area  available  for 
producing  new  ring  fractures.  Although  overlapping  ring  frac¬ 
tures  do  occur,  it  is  assumed  here  that  once  a  ring  fracture 
is  formed  there  is  no  possibility  for  another  ring  to  be 
formed  in  the  same  area.  This  is  a  valid  assumption  for  the 
shallow  fracture  surfaces  where  an  impact  in  the  vicinity  of 
a  pre-existing  ring  is  more  likely  to  contribute  to  the  growth 
of  an  erosion  pit  than  the  formation  of  a  new  ring  fracture. 

It  is  for  this  reason  that  the  total  effective  area  of  the 
ring  fractures  present  at  any  instant  is  also  excluded  from 
the  surface  area  available  for  the  formation  of  ring  fractures. 

Referring  to  Fig.  4  and  5  it  can  be  seen  that 
very  few  ring  fractures  are  produced  during  the  earliest  stage 
of  the  erosion  process.  There  is  abundant  evidence  that 
numerous  particles  have  struck  the  surface  by  the  distribu¬ 
tion  of  very  small  nicks  in  the  surface,  however,  these  flaws 
exhibit  relatively  minor  growth  during  subsequent  exposure 
to  the  particle  impacts.  The  ring  fracture  population  begins 
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to  Increase  after  an  incubation  period.  It  is  conjectured 
that  ring  fractures  are  being  produced  on  the  surface  due  to 
multiple  impacts  at  velocities  which  are  substantially  below 
those  which  would  be  required  to  produce  a  ring  fracture  due 
to  a  single  particle  impact.  The  reasoning  is  as  follows. 

The  original  surface  has  an  initial  flaw  distribution,  but  as 
the  particles  collide  with  the  surface  a  more  extensive  flaw 
distribution  evolves.  The  probability  that  a  ring  fracture 
will  result  is  increased  with  increasing  exposure  time.  This 
effect  can  be  seen  most  clearly  in  Fig.  5  where  a  negligible 
number  of  ring  fractures  are  observed  on  the  surface  of  the 
specimen  after  an  exposure  time  of  one  minute.  During  the 
next,  one  minute  increment  in  the  exposure  time  numerous  ring 
fractures  are  formed.  Research  on  the  evaluation  of  the  pro¬ 
bability  for  ring  fracture  formation  on  glass  plates  due  to 
solid  particle  impacts  has  concentrated  on  single  particle 
collisions  over  a  range  of  particle  sizes  which  is  above  that 
usually  of  interest  in  erosion  studies.  In  order  to  account 
for  the  observed  time-dependent  probability  function  in  multiple 
particle  collisions,  it  is  necessary  to  devise  an  adjunct 
experimental  program  to  obtain  the  data  required.  Considera¬ 
tion  is  being  given  to  static  and  dynamic  indentation  tests 
using  particle  sizes  corresponding  to  those  used  in  the  erosion 
tests  and  glass. specimens  which  have  been  exposed  to  multiple 
particle  impacts  for  a  range  of  exposure  times  which  are  less 

than  those  at  which  ring  fractures  are  first  observed.  It  is 
(?) 

knownv  '  that  the  threshold  fracture  stress  in  dynamic  tests 
is  higher  than  that  found  in  static  tests,  however  the  attain¬ 
ment  of  the  velocities  required  for  single  impacts  of  glass 
beads  less  than  300  microns  in  diameter  would  involve  a  much 
more  extensive  experimental  facility  than  the  static  tests. 

Let  (p  be  the  probability  that  a  ring  frac¬ 
ture  will  be  produced  under  conditions  correspond ing  to  those 
for  a  lightly  eroded  surface.  The  time  lag  in  the  formation 
of  ring  fractures  on  the  surface  of  the  glass  specimens  can  be 
represented  by 
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with  the  values  of  the  parameters  m][  and  mg  being  obtained 
from  microscopic  observations  on  ring  fractures.  General 
trends  in  <p(t)  for  different  size  beads  and  impact  velocities 
can  be  noted  in  order  to  establish  relationships  between  m^, 
trip  and  the  variables  which  are  to  be  specified  in  the  model. 

As  material  begins  to  be  removed  from  the 
surface  layer,  an  increasing  number  of  impacting  beads  con¬ 
tribute  to  the  growth  of  the  nucleated  erosion  pits  and  cannot 
be  considered  to  be  available  for  ring  fracture  production. 

In  order  to  take  this  effect  into  account  in  the  modeling 
studies,  the  form  of  the  time-dependent  ring  fracture  proba¬ 
bility  function  in  Eq.  (6)  can  be  modified  so  that  after 
reaching  a  maximum  it  approaches  zero  as  the  exposure  time 
increases.  It  is  preferable  to  modify  <p  and  let  the  number 
of  impacting  particles  remain  constant,  instead  of  trying  to 
devise  a  method  for  deciding  on  the  number  of  beads  contribu¬ 
ting  to  ring  fracture  formation  and  the  number  contributing 
to  the  growth  of  erosion  pits  as  a  function  of  time. 

The  formation  of  erosion  pits  is  a  function 
of  the  ring  fractures  formed  on  the  surface  of  the  specimen. 
Guided  by  the  microscopic  investigations  of  multiple  particle 
impacts, 

Nu  =  g(D)N(t)  (7) 

where  Nu(t)  is  the  number  of  erosion  pits  nucleated  at  time 
t  per  unit  surface  area 

g(D)  is  a  function  of  the  thickness  of  the  surface 

layer  which  is  determined  by  the  penetration  of 
the  fracture  surface  into  the  bulk  of  the  speci¬ 
men  and  is  dependent  on  the  radius  and  velocity 
of  the  impacting  glass  beads. 
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For  the  very  small. 
Each  r inp  fracture 
»  r<  si  on  pit  a;;  can 


glass  heads,  g(D)  can  be  equated  to  unity, 
j  fr.sent  ially  a  nucleation  site  for  an 
to  seen  in  'Kip.  4  and  5. 


The  area  removed  from  the  surface  of  the 
specimen  through  the  growth  of  erosion  pits  is  specified  by 

f  t  d  N  ( r  ) 

aAl(t)  =/o  ~dT -  Ki(t-r)dT  (8) 

which  simply  is  the  summation  of  the  number  of  erosion  pits 
nucleated  in  the  t ime  interval  r,  T+dr  times  the  growth  rate, 
K,,  they  experienced  fiom  the  time  they  were  formed  to  the 
current  time  t. 


Combining  Eq.  (5),  (7),  and  (8), 

=  J  </>(t  )dT-g(D)J  <p(r)  I  d-^-Le-I  K1(T-e)dedT-Ao/  N(T)<p(T)dr  (9) 

The  erosion  of  the  surface  layer  can  be  evaluated  by  solving 
this  single  integral  equation  for  N(t).  The  erosion  tests 
and  adjunct  experimental  programs  are  being  used  in  conjunc¬ 
tion  with  the  Hertzian  analysis  of  solid  particle  impacts  to 
supply  explicit  representation  for  the  functions  <p(t),  K^(t), 
g(D),  and  AQ. 

A  limiting  case  of  Eq.  (9)  can  be  considered 
to  substantiate  its  validity.  For  very  small  beads  at  a  suf¬ 
ficiently  high  impact  velocity  the  probability  for  ring  frac¬ 
ture  formation  approaches  unity  and  g(D)  should  also  he  unity, 
then  Eq.  (9)  reduces  to 

-  f  d-r  -  f  f  K1fr-e)  dedT-AQ  f  N(r)dr  (10) 

-'O  b  •'O  'O 

To  further  simplify  this  equation,  assume  K^(t.)  =  0  for 
'■)*rt<oo,  which  means  -  0,  then 
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■ggr  +  nA0N  =  n 

N(t)  =  (11) 

The  maximum  number  of  ring  fractures  which  can  be  found  on  a 
unit  area  of  the  surface  of  the  first  platelet  is  given  by 


N 

max 


(12) 


where  Ac  is  the  area  of  a  single  ring  fracture.  It  is  seen 
that  this  is  a  reasonable  result  in  view  of  the  assumptions 
made.  The  form  of  Eq.  (11)  also  shows  that  the  number  of 
ring  fractures  increases  as  a  function  of  time  when  the  number 
of  impacts  on  the  surface  is  increased. 


Once  AA-^(t)  is  determined  from  Eq.  (8),  it 
can  be  substituted  into  Eq.  (4)  to  obtain  the  volume  eroded 
from  the  first  platelet.  It  was  stated  earlier  that  the 
eroded  volume  from  the  first  platelet  was  composed  of  two  com¬ 
ponents:  AV1  and  AV^.  As  AV^  increases  with,  time,  a  distri¬ 
bution  of  conical  frustrums  are  exposed  on  the  eroded  surface-. 
This  phenomenon  is  particularly  evident  for  the  290  micron  beads 
as  shown  in  Fig.  6(b).  The  frustrum  density  after  removal 
of  the  interstitial  material  reaches  a  fairly  saturated  state 
before  the  frustrums  are  eroded  further.  The  distribution  of 
frustrums  on  the  specimen  surface  constitute  the  portion  of 
the  first  layer  for  which  AV^  represents  the  eroded  volume. 

The  parameter  rj  is  taken  to  be  a  constant  defined  by 


ii»AVi(t)=^1  (13) 

Then  the  evaluation  of  tj  can  be  mqde  by  assuming  that  a  typi¬ 
cal  frustrum  is  enclosed  by  a  rectangular  parallelepiped  with 
2 

volume  DH  .  The  parameter  r?  is  related  to  the  difference  in 
volume  between  the  parallelepiped  and  the  volume  of  the  frus¬ 
trum, 
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where  r  is  the  radius  of  the  ring  fracture  at  the  surface 
°  2 

of  the  platelet  (A0=7rro  )  and  S  is  the  radius  at  the  base  of 
the  frustrum  a  depth  I)  below  the  surface.  It  is  observed 
that  a  near  dense  array  of  frustrums  remains  after  the  inter¬ 
stitial  material  is  removed  for  the  290  micron  beads.  In  this 
case  H»2S,  so 

r  r  2 

V  ~  1  -  Y2  l+~  +  (~g~)  (15) 

Microscopic  observation  of  the  erosion  pro¬ 
cess  as  a  function  of  time  for  the  290  micron  beads  shows 
that  the  density  of  ring  fractures  increases  to  a  maximum 
value  and  then  decreases  toward  a  number  approaching  that  for 
an  array  of  frustrums  whose  bases  almost  completely  cover  a 
unit  area.  Frustrums  are  removed  with  the  interstitial  mate¬ 
rial  when  their  bases  are  not  in  contact  with  the  bulk  mate¬ 
rial  but  lie  in-between  two  adjacent  frustrums  ^ .  This 
behavior  should  be  taken  into  account  in  the  present  formula¬ 
tion  of  the  model.  In  the  case  of  the  70  micron  beads,  the 
growth  of  nucleated  pits  and  chipping  of  conical  frustrums 
occur  a  very  short  time  after  a  ring  fracture  develops. 

The  form  of  the  governing  equation  for  the 
erosion  of  the  surface  layer  in  the  basic  model  can  be  com¬ 
pared  with  the  corresponding  equation  derived  in  the  earlier 
stages  of  this  research  (Reference  1,  Eq.  (4))  which  may  be 
written 

AV1(t)  =  D  j  G^(t-T)2  I1(r)(l-AA1(T))d r  (l6) 

Substituting  Eq.  (7)  into  (8)  and  using  Eq.  (5),  Eq.  (4) 
becomes 

AV1(t)  =  r?Dj  K1(t-T)g(D)nv’(r)  (l-AA1(r  )-N(r)A0)dr  (17) 
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p  p 

The  general  growth  rate  function  K^(t-r)  replaces  G^  ft-r) 
and  the  pit  nucleation  probability  function  is  now  more 
explicit, 

^(t)  =  g(D)n*(t)  (18) 

The  microscopic  observation  of  the  eroded  surfaces  also  showed 
that  the  area  of  the  number  of  rings  nucleated  should  be 
excluded  from  the  area  available  for  additional  ring  fracture 
formation  as  Indicated  by  the  last  term  in  Eq.  (17). 

Now  that  the  erosion  mechanisms  contributing 
to  Av?  and  AV^  have  also  been  identified,  more  explicit  mathe¬ 
matical  representations  will  be  provided  for  determining  the 
volume  removal  due  to  these  stages  of  the  general  erosion  pro¬ 
cess.  For  different  erosion  conditions  the  measured  differences 
in  the  weight-loss  data  are  associated  with  the  time  at  which 
each  erosion  mechanism  (contributing  to  AV^,  AV^,  or  AV^) 
becomes  operative  and  the  degree  to  which  it  influences  the 
overall  erosion  behavior.  It  is  now  concluded  that  the  over¬ 
all  erosion  behavior  can  be  described  for  a  broad  range  of 
bead  diameters  and  impact  velocities  in  terms  of  the  present 
formulation  of  the  analytical  model. 

3.2.3  Experimental  Evaluation  of  Model  Parameters 

Programs  are  underway  to  acquire  data  per¬ 
taining  to  each  of  the  parameters  required  to  obtain  explicit 
results  from  the  analytical  model  as  represented  by  Eq.  (9), 

Each  of  the  four  parameters,  A0,  <p,  D,  and  will  be  dis¬ 
cussed  here  In  relation  to  methods  for  obtaining  the  necessary 
infer  mat  ion. 


3. 2, 3.1  Ring  Fracture  Si3e  Distribution  and 

Probability  Function 

Quantitative  data  pertaining  to  the  ring 
fractures  formed  by  the  collision  of  the  glass  beads  with 
glass  plates  on  a  scale  appropriate  to  solid  particle  erosion 
are  required  in  the  development  of  the  basic  model.  Tests 
PY-12  to  PY-19  (see  Table  I)  were  run  to  obtain  information 
on  the  probabilities  for  ring  fracture  formation  as  a  function 
of  impact  velocity.  In  order  to  determine  the  threshold  velo- 
cit3'-  for  ring  fracture  initiation,  specimen  PY-12  was  run  for 
10  sec.  intervals  beginning  at  38  fps  and  proceeding  in.  incre¬ 
ments  of  19  fps  in  the  velocity  until  ring  fractures  were 
observed  on  the  surface.  When  ring  fractures  were  found  at 
approximately  170  fps,  the  number  of  rings  at  several  loca¬ 
tions  on  the  specimen  surface  were  counted  using  an  optical 
microscope  at  440X.  These  areas  were  then  scanned  a  second 
time  and  the  diameters  of  the  rings  were  measured  and  recorded. 
The  distribution  of  ring  diameters  for  impact  velocities  of 
,188,  208,  226  and  246  fps  are  given 'in  Fig.  7  to  10.  The  - 
avera ge  diameter  in  each  case  is  listed  in  Table  VI.  Ihe 
average  diameter  of  the  ring  fractures  for  an  impact  velocity 
of  246  fps  is  decisivply  greater  than  226  fps  and  the  distri¬ 
bution  curves  show  a  shift  to  an  increased  population  of 
larger  diameter  rings.  Each  value  was  determined  on  a  samp¬ 
ling  of  200  rings.  The  agreement  with  the  Hertzian  theory 
(Table  III)  is  excellent.  Comparing  the  average  ring  diameter 
for  lyrex  impacted  by  the  290  micron  beads  at  approximately 
200  fps  with  the  corresponding  value  for  fused  silica,  3. 17 . 
mils  compared  with  2.63  mils^,  it  is  seen  that  there  is  a 
sizable  difference  between  the  two  glasses.  The  explanation 
for  this  difference  requires  further  investigation. 

h  • .  The  near  equality  of  the  average  ring  dia¬ 

meter,  at  velocities  below  226  fps  is  not  given  any  signifi¬ 
cance  at  this  time,  since  the  sampling  at  208  fps  was  130 


Figure  7.  Distrib 
for  290-m: 


2  31-2.56  2.56-2.81  2.81-3.06  3.063.31  3.31-3.56  3.563.81  3.81-4.06  4.064.31  4.31-4.56 

Ring  Fracture  Diameter  (mils) 

FIGURE  8.  Distribution  of  ring  fracture  dia.neters 
for  290-micron  beads  at  208  fps 
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’racture  diameters 
246  fps 


TABLE  VI 

RING  FRACTURE  PROBABILITIES  AND  AVERAGE  DIAMETERS 


Impact 

Velocity 

Number  of 
Impacts 
per  sq.  in. 

Number  of 
Ring  Fractures 
per  sq.  in. 

Average 
Diameter 
of  Ring 
Fracture 
(mils ) 

RPM 

fps 

200 

188 

9,500 

1,725 

0.18 

3.15 

220 

208 

10,400 

5,710 

0.55 

3.17 

240 

226 

11,400 

7,900 

O.69 

3.17 

260 

246 

12,  300 

11,000 

0.90 

3.38 

rings  and  at  18b  fps  it  was  only  40  rings  due  to  the  sparse 
ring  fracture  population.  These  tests  will  be  continued  at 
velocities  greater  than  those  reported  here  to  observe  further 
the  trend  in  the  average  ring  diameter  as  a  function  of  velo¬ 
city.  So  far  the  results  are  most  promising  in  terms  of  the 
Hertzian  analysis,  since  ring  fractures  generally  do  occur 
slightly  beyond  the  contact  area  as  predicted  by  Hertz  theory. 

It  was  already  pointed  out  that  the  proba¬ 
bility  for  the  formation  of  ring  fractures  is  a  time-dependent 
function  due  to  the  changing  surface  condition  of  the  glass 
plate  when  progressively  bombarded  by  the  glass  beads.  Quan¬ 
titative  measurements  of  the  t ime-variation  of  the  probability 
function  <p  will  be  made  subsequently.  The  tabulations  in 
Table  VI  are  simply  based  on  the  number  of  ring  fractures 
produced  by  a  known  number  of  bead  impacts.  These  results 
are  still  tentative  but  do  provide  some  perspective  on  the 
relationship  between  the  impact  velocity  and  9. 

A  series  of  tests  (PY-20  and  PY-21)  were 
run  in  order  to  obtain  the  distribution  of  ring  fracture  dia¬ 
meters  for  the  70-micron  beads  (Fig.  11  and  12),  The  mean 
diameter  at  375  fps  is  1.25  mil  and  1.28  mil  at  4-50  fps.  The 
mean  depth  of  the  ring  fracture  surface  was  found  to  be 
approximately  9  microns  at  375  fps  and  10  microns  at  450  fps. 
Both  values  are  consistent-  wi»h  the  Hertzian  predictions  of 
the  maximum  contact  diameters  which  are  slight ly  less  than 
1  mil  as  computed  in  Table  ITT.  The  data  obtained  are  based 
on  the  measurement  of’  200  ring  fractures  using  optical  micro¬ 
scopy  at  a  magnification  of  4407. 


3. 2. 3. 2  Ooong’fry  '  r  ’ 'on  Ira  I  t  car  <  u  On  r  fares 

V/hr-Ti  v  t  ewl  fit  t  hr  rri.drd  nurf'nrr  of  t  hr 
glass  specimens  in  t  hr  ml  mxwi.pr ,  '  t , .  Irp'  f  "t  1  hr  eoh 
fractures  can  tie  measured  In  ml  rmnn  1  1  ■ 1 ' ,i!'  a 1 
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0.86-0.98  0.98-1.10  1.10-1.22  1.22-1.34  1.34-1.46  1.46-1.58  >1.58 

Ring  Fracture  Diameter  (mils) 

FIGURE  Distribution  of  frng  fracture  diameters  for  70-micron  glass 
beads  impacting  or.  glass  at  450  fps. 


surface  and  then  the  base  of  the  frustrums  once  they  are 
exposed.  The  lateral  scale  on  the  microscope  can  be  read  to 
an  accuracy  of  a  few  microns.  The  tangent  of  the  cone  semi¬ 
angles  obtained  in  this  way  never  agreed  with  the  values  of 
68.  5°  for  glass  and  05. 5°  for  fused  silica^^  measured  under 
quasistatic  loadings.  The  measured  values  of  the  semiangles 
under  the  present  impact  conditions  were  found  to  range  from 
45°  to  60° .  A  similar  observation  was  made  by  Bowden  and 
Field'  '  in  conjunction  with  a  lead  slug  striking  the  edge 
of  a  0.25  in.  glass  plate  which  was  3  inches  square.  They 
observe  from  their  high-speed  photographic  studies  that  an 
important  feature  of  cone  fractures  is  their  angle  of  inclin¬ 
ation  with  the  top  surface  of  the  glass  plate  which  lies 
primarily  in  the  range  of  50  to  60  degrees.  This  means  that 
the  included  semi-angle  is  in  the  range  of  30  to  40  degrees 
which  is  in  conformity  with  the  rough  estimates  made  for  the 
semi-angles  for  the  primary  cone  fractures  investigated  in 
this  program.  Bowden  and  Field  note  that  if  the  horizontal 
and  vertical  displacements  at  the  surface  for  the  steady  state 
solution  of  a  Rayleigh  wave  are  used  to  compute  the  direction 
of  the  maximum  tensile  stress,  an  included  semi-angle  of  3^° 
is  predicted.  Due  to  the  passage  of  a  Rayleigh  wave  the  mate¬ 
rial  particles  at  the  surface  of  a  semi- infinite  half-space 
follow  an  elliptical  path  with  major  axis  in  the  vertical 
direction.  When  Poisson's  ratio  is  0.25,  it  turns  out  that 
the.  cotangent  of  the  ratio  of  the  major  to  minor  axis  of  the 
elliptical  path  corresponds  to  an  angle  of  approximately  34°; 
however  the  pnysical  meaning  of  this  result  is  not  clear. 

A  number  of  investigators  have  considered 
the  formation  of  cone  fractures  in  glass  plates,  however  there 
are  only  limited  data  pertaining  to  the  fracture  of  glass 
specimens  impacted  by  small  diameter  solid  spheres.  While 
qualitative  aspects  of  the  fracture  process  under  these  con¬ 
ditions  can  be  described,  the  quantitative  evaluation  of  the 
derived  relations  must  await  future  experimentation  on  the 
specific  materials  used  in  the  erosion  tests. 
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The  experimental  work  of  Hamilton  and 
Rawsonw/  on  two  different  glasses  with  various  surface  treat¬ 
ments  indicates  that  the  linear  proportionality  betvreen  the 
applied  load  and  the  indentor  radius  at  fracture  is  not  gen¬ 
erally  true  in  the  Auerbach  range  The  following  relation 

is  obtained 

Pc  =  KRm '  (19) 

where  the  values  of  K  and  m1  listed  in  Table  VII  are  deter¬ 
mined  from  the  data  plotted  in  their  paper.  PQ  is  the  cri¬ 
tical  fracture  load  applied  to  a  spherical  indenter  of  radius 
R.  The  values  of  the  parameters  K  and  m'  are  for  Pc  in  kilo¬ 
grams  and  R  in  millimeters.  On  the  basis  of  the  variability 
in  the  form  of  Eq.  (19),  Hamilton  and  Rawson  argue  that  the 
linear  P  -R  (Auerbach)  relation  is  not  a  universal  relation 
and  thereby  dismiss  Roesler's^  argument  that  the  linearity 
of  the  P  -R  relation  would  require  a  special  flaw  distribution 
function  common  to  all  brittle  surfaces.  It  is  our  observa¬ 
tion  from  the  data  reported  by  Hamilton  and  Rawson  that 
the  relation  for  the  Hertzian  range 

P  =  k’r2  (20) 

C 

is  similarly  replaced  by  a  general  relation  of  the  form  given 
in  Eq.  (19).  The  values  of  K  and  m'  for  this  case  are  also 
listed  in  Table  VII. 

Roesler^^  developed  a  scaling  relation 
for  describing  the  propagation  of  a  Hertzian  crack  as  the  load 
applied  to  a  cylindrical  indenter  is  increased.  The  radius 
of  the  ring  fracture,  S,  measured  parallel  to  the  plane  of 
the  glass  plate  asymptotically  reaches  a  state  which  is  des¬ 
cribed  by  the  relation 


2S 


Kp2/3 


(21) 
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TABLE  VII 

EXPERIMENTALLY- EVALUATED  PARAMETER  FOR 
THE  DETERMINATION  OF  THE  CRITICAL  FRACTURE  LOAD 


Auerbach 

Range 

Hertz 

Range 

m' 

K 

m' 

K 

Polished  Plate  Glass 

As  received 

1.032 

14. 1 

1.706 

2.28 

4  minute  etch 

0.812 

32.0 

1.671 

5.29 

7  minute  etch 

0.720 

43.6 

1.558 

36.2 

Float  Glass 

Tin  bath  surface 

0.778 

35.5 

1.606 

6.14 

Free  surface 

0.694 

43.0 

1.576 

37.7 

Roesler  estimates  that  this  relation  holds  for  a  load  P  equal 

to  twice  the  load  required  for  initiation  of  fracture,  P  , 

c* 

since  after  this  point  crushing  of  the  glass  surface  takes 
place.  The  experiments  in  support  of  the  scaling  relation 
derived  by  Roesler  were  carried  out  with  a  flat-faced  cylin¬ 
drical  indentor  instead  of  a  spherical  indentor.  According 
to  Roesler  these  test  conditions  were  used  rather  than  the 
customary  spherical  indentor,  since  spherical  indentors  do 
not  produce  large  and  regular  cone  fractures.  For  spherical 
indentors  secondary  fractures,  multiple  ring  fractures,  de¬ 
velop  as  the  contact  area  grows  with  increasing  load.  It 
will  be  assumed  that  Eq. (21)  is  also  valid  for  a  sphere 
impacting  an  elastic  half-space. 

Benbow^)  carried  out  a  series  of  experi¬ 
ments  employing  the  cylindrical  indentor  and  an  unidentified 
fused  silica.  The  experimental  results  are  reported  in  gra¬ 
phical  form  in  his  paper.  The  values  obtained  from  the  line 
fitted  to  the  experimental  data  provides  an  explicit  repre¬ 
sentation  for  Eq.(21).  It  is  estimated  that  k  =  3.84  x  10“  3 
in/(lb)2/3. 


The  depth  of  penetration  of  a  single  cone 
fracture  can  be  determined  from  the  information  which  has 
been  accumulated .  The  load  applied  to  the  specimen  by  an 
impacting  bead  is  determined  by  equating  the  kinetic  energy 
of  the  impacting  particle  to  the  strain  energy  at  the  time  of 
maximum  impression.  According  to  the  Hertz-Huber  theory, 


^  /5 


3/5 


Pm  "  3  +  k2^ 


-2/5  6/5  2 
v_  R 


(22) 


where  p2  is  the  density  of  the  impacting  sphere 

V  is  the  impact  velocity 
°  p  p 

^  _  — =-i— ,  k0  =  '  —  r-f ■  * ~  are  the  elastic  moduli  for  the 

1  E!  d  b,2 

sphere  and  half-space,  respectively,  (v  is  Poisson's 
ration  and  E  is  Young's  modulus). 


Values  of  %  are  listed  in  Table  III  for  the  materials  and 
impact  conditions  used  in  the  experimental  program.  Assum¬ 
ing  Bossier' s  scaling  relation  is  valid  for  a  spherical 
indenter. 


2S  =  i(P 


when  P_  >  P„ 
m  c 


where  according  to  the  experimental  data  of  Hamilton  and 
Raws  on  the  critical  fracture  load  is  given  by  Eq. (19). 
Hence,  referring  to  Fig.  13»  we  finally  obtain 


tan/3 


(24) 


where  r  =  K„a  is  the  radius  of  the  inner  ring  fracture  at 
o  2 

the  surface  of  the  specimen  ana  a  is  the  radius 
of  the  contact  area  computed  from  Herts  theory.. 


For  the  present  only  very  rough  approxi¬ 
mations  will  be  made  for  the  parameter  appearing  in  Eq.(24). 
The  coefficient  k  in  Eq.(23)  will  be  equated  to  the  value ..... 
determined  by  Benbow  for  fused,  silica:  k  =  3.84  x  lCT%n/  -■ 
(lb)2//3.  The  radius  of  the  inner  fracture  ring  will  b.e 
equated  to  the  radius  of  the  maximum  contact  area  deter¬ 
mined.  from  the  Hertz-  H-’ber  theory  and  listed  in  Table  III; 
and,  until  more  accurate  measurements  have  been  made,  /3  will 
be  set  equal  to  60° .  The  values  of  D  calculated  for  both 
bead  radii  and  .listed  in  Table  VIII  over  the  velocity  range 
used  in  the  experimental  program  are  in  general  agreement 
with  the  experimentally-determined  values,  Table  VIII  also 
provides  the  values  of  the  shape  factors,  rj ,  calculated  from 
Eq.  (15)  using  the  computed  value  of. ..S  and  the  value  of  rQ=a 
as  recorded  in  Table  III. 


•  ■-  TABLE-  VIII 

COMPUTED  FRACTURE  CONE  GEOMETRIC  PARAMETERS 


Radius  at 

Depth  of  cone 

"•  Shape 

base'  of  cone 

frabture  below 

factor 

v 

fracture,  S ( in . ) 

surface,  D  (in, ) 

n 

70  Micron _ 
Glass  Beads 

• 

0.204  x  10"3(  5.2#*) 

200  fps 

0.73  y  10” 3 

0.552-; 

300  fps  v  . 

0.96  X  10’ 3 

0.306  x  10’3(  7.8m) 

0.577 

400  Ops 

1.22  X  10" 3 

0.433  x  10" 3 (  11.0m) 

0.601 

500  fps 

1.44  X  10" 3 ..." 

0.508  x  10” 3 (  12.9M) 

0.614 

290  Micron 
Glass  Beads 


200  fps 

4.48  x 

CO 

1 

G 

rri 

1.75 

x  10" 3  (  44,5  m) 

0.628 

300  fps 

-  6,25  x 

10" 3 

2.62 

x  10"  3  (  66.5m) 

0.648 

400  fps 

7.90  X 

10‘3 

3.46 

x  TO* 3 (  88.0  m) 

0.660 

500  fps 

9.35  x 

10"3 . . 

4.20 

x  10" 3  (106.5  m) 

O.667 
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;  '  3.2. 3*3  "F^.oded  Area  Growth  Rates- 

,xx  The  growth  rate,  K-jJt),  for  nucleated 
erosion  pits  can  be  determined  through  direct  observation  of 
areas  from  which  material  is  being  removed  as  a  function  of 
the  exposure  | ime ;  however,  for  both  bead  sizes  the  areafe.  ■ 
.patterns? have  geometrically  intricate  forms  which  are  tedious 
to  measure  accurately.  In  Order  to  obtain  more  detailed  : 
information  for  the;  Statistical  analysis  of : the  growtiffrate  % 
function,  an  electronic"  imag§  "  analyzer  was  used:  in.  conjunc¬ 
tion  with  a  low  power  optical  mi'cr osbbpe .  ' Tfte'lsxper  imental 
setup  was  only  preliminary  in  order  to  evaluate  the'  feasi-  ' 
bility  of  using  this  arrangement  to  obtain  quantitative  mea¬ 
sures  of  the  change  in  the  eroded  area  at  a  number  of  locations 
on  the  specimen  surface  as  a  function  of  the  time  of  exposure. 


,The  evolution  of  damage  on  the  surf ce  of  a 
glass  specimen  impacted- by  70  micron  beads  is  shown,  in  Fig.-  ; 

4  and  5.  In  bright  field. the  details  of  the  eroded  surface 
are  suppressed  but  the  outline  of  regions  from  which  material 
has  been  removed  standout  as  black  patches  on  a  white  surface. 
A  video  camera  is  used  to  view  a  portion  of  the  specimen  sur¬ 
face  through  the  eyepiece  of  an  optical  microscope,  A  magni¬ 
fication  factor  of  15  is  satisfactory  in  accordance  with  the 
scale  of  the  erosion  damage.' 

-  By  using. the  specular  reflection  charac¬ 
teristics  of  a  smooth  reflective  surface,  measurements  can 
be  made'  of  the  percentage  of  an  arbitrary  rectangular  area 
that  is,  undamaged.  The  first  step  in  this  evaluation  is  to 
illuminate  the  surface  to  separate,  through  the  reflected 
light  intensity,  the  original  surface  areas  from  the  'pitted 
arbas.  This  is  done  by  forcing  the  illumination  angle  and 
the  viewing  angle  to  both  be  perpendicular  to  the  surface. 

The  undamaged  areas  then  appear  very  bright  due  to  the  strong 
specular  reflection  from  them,  while  pitted  areas  appear  dark. 
This  lighting  scheme  is  depicted  in  Fig.  14.  If  the  intensity 


Figure  it.  Illumination  scheme  to  separate  damaged  regions  from 
undamaged  regions  by  intensity 


Relative 

Area 


Figure  15  Intensity  distributions  of  partially  damaged  surface 
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'•  Li. 


distribution  of . the  surface,  when  viewed  with  this  lighting 
were  graphed,  it  would  appear  as  in  Pig.  15. 


In.  subsequent  development  of  the  areal 


measurement^  equipment  the  surface,  illuminated  in -the  manner 
'described,  will  be  viewed  by  the  television  camera  with  the 
aid  of  a  microscope  equipped  with  a  two-dimensional  micrometer 
^ stage,  so  the  same  area  can  be  vlswed  after  additional  expo- 
to  the  erosive  environment.  Since  the  video  signal  is 
^  a {faithful  voltage- analog  of  the  light,,  intensity,  the,  voltage 


distribution  of  the  signal  is  the  same  as  that  shown  in  Fig. 


14. 


This,  video  signal  is  first  split  into  two 
components.  A  threshold  circuit  has  a  bilevel  output;  it.  is 
high  whenever  the  input  signal  Is  above  an  arbitrary  but  fixed 
level  and  low  when  the. input  signal  is  below  that  level.  If 
this  threshold  level  is  chosen  to  be  T  as  shown  in  Fig.  15, 
then  the  threshold  output  will  be  high  when  undamaged  surface 
Is  being  scanned  by  the  television  camera  and  low  when  the 
surface  being  scanned  is  damaged.  The  relative  amount  of  time 
the  signal  is  high  is  then  proportional  to  the  amount  of 
undamaged  surface  area  in  the  video  picture.  The  rest  of  the 
processor  consists  of  a  circuit  to  measure  the  amount  of  time 
the  threshold  signal  is  high,.  This  is  accomplished  by‘ gating 
an  oscillator  to  a  counter  with  the  threshold  video.  An 
additional  useful  feature  is  a  variable  window  which  can  b  e 
used  to  frame  exactly  that  region  to  be  measured. 

The  arrangement  used  in  the  initial  tests 
included  a  plate  which  masked  the  face  of  the  specimen  except 
for  a  hole  approximately  0.2  inch  in  diameter.  A  holder  was 
fabricated  for  the  specimen  and  a  scale  was  used  to  position 
the  hole  over  the  same  area  bn  the  specimen  at  the  next  incre¬ 
ment  in  the  exposure  time.  The  accuracy  in  returning  to  the 
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same  location  was  not  satisfactory,  so  the  readings  taken  at 
this  stage  in  the  development  of  the  equipment  are  of  little 
value.  The  sides  of  the  hole  in  the  plate  also  produced  a 
reflection  on  the  video  picture  which  further  hindered  the 
attainment  of  acceptable  results. 


i.  • ' 


3.2.4  'Weight-loss.  .Data 

Weight-loss  data  were  obtained  for  the  direct 


impact  of  the,. 290  micron  glass  beads  in  test  numbers  •PY^Q  and 
PY-11,  and  for  the  70  micron  glass  beads  in  test  numbers  PY-22 
and  PY-23  as  recorded  in  Table  I.  The  weight- loss  curves  in 
each  case  are  shown  in  Fig.  16  to  19.  Reference  to  Table  V 
provides  information  on  the  area  exposed  to  the  erosive  envi¬ 
ronment  and  estimates  of  the  number  of  particles  striking  the 
exposed  area.  It  should  be  noted  that  the  direct  impact  tests 
for  the  290  micron  glass  beads  were  for  a  different  experi¬ 
mental  arrangement  than  that  shown  in  Fig,  1.  In  this  case 
2.25  sq.in.  of  the  specimen's  surface  are  exposed  to  the  bead 
impacts.  In  these  initial  tests,  the  arrangement  in  the  AFML-  . 
Bell  erosion  facility  for  dropping  glass  beads  in  the  path  of 
the  specimen  at  the  end  of  the  rotating  arm  resulted  in  a  non- 
uniform  distribution  of  beads  across  a  portion  of  the  lateral 
dimension  of  the  specimen.  The  disadvantage  of  this  arrange¬ 
ment  is  that  different  areas  on  the  specimen  surface  will  be 
at  different  stages  of  the  erosion  process  duo  to  the  non- 
uniform  number  of  bead  impacts  occurring  along  the  specimen 
surface.  The  measured  weight- Toss  data  are  then  reflecting 
the  cumulative  effect  of  two  or  more  erosion  mechanisms  being 
operative  as  a  function  of  the  lateral  dimensions  of  the 
specimen  and  not  simply  a  function  of  the  exposure  time-.- 
Simple  modifications  were  made  in  the  nozzle  at  the  end  of  V 
the  blast  tube,  but  the  desired  result  was  not  achieved.. 

Cover  plates  were  th< n  used  as  indicated  in  Fig.  1  so  that  the 
specimen  is  exposed  to  a  uniform  distribution  of  particles 
which  exists  near  the  center  of  the  stream  of  particles 
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falling  in  the  path  of  the  specimen.  This  modification  in 
the  testing  procedure  is  more  in  line  with  the  analytical 
model  which  does  account  for  a  temporal  but  not  a  spatial 
dependence  of  the  erosion  mechanisms  that  have  been  identified 
in  the  solid  particle  erosion  of  glass  specimens.  The  nozzle 
used  in  tests  PY-10  and  PY-11  did  spread  the  beads  over  the 
full  length  of  the  specimens;  but  when  the  arm  was  rotating, 
the  density  of  bead  collisions  with  the  specimen  was  skewed 
to  the  outboard  end. 

The  general  form  of  the  weight- loss  curve 
in  Fig.  16  correlates  quite  well  with  the  microscopic  obser¬ 
vations  of  the  operative  erosion  mechanisms  as  described  in 
Section  3.2.1,  On  the  basis  of  the  analysis  provided  in 
Section  3. 2. 3*2  it  is  also  possible  to  check  quantitative 
details  of  the  erosion  process.  The  cone  depth  calculated 
for  290  micron  beads  at  an  impact  velocity  of  200  fps  (test 
PY-10)  is  1.75  mils  (Table  VIII).  The  mass  of  material  com¬ 
prising  the  surface  layer  of  the  specimen  as  defined  in 
Section  3.2,2  is  144  mg.  The  volume  of  the  surface  layer  is 
divided  into  two  components:  AV^  representing  the  interstitial 
material  for  a  dense  rectangular  array  of  frustrums  covering 
the  exposed  surface  area  of  the  specimen,  and  AVg  corres¬ 
ponding  to  the  volume  contained  within  the  fru strums.  The 
mass  of  the  material  in  AV^  is  90  mg.  using  the  shape  factor 
in  Table  VIII.  The  curve  in  Pig,  16  shows  an  initial  incuba- . 
tion  period  before  measurable  weight  loss  occurs,  a  period  of 
accelerated  weight  loss  when  the  material  is  removed  from  the 
highly-fractured  interstitial  region  between  the  more  erosion, 
resistant  frustrums,  and  a  period  of  uniform  material 
removal  corresponding  to  chipping  of  the  frustrums  and  the 
bulk  material.  Microscopic  examination  of  specimen  PY-10  at 
the  end  of  the  test  indicated  that  most  of  the  surface  layer 
was  removed  except  for  a^sma 11  portion  of  the  surface  at;  the 
inboard  end  of  the  specimen  where  frustrums  were  still  intact. 
The  calculated  mass  of  material  corresponding  to  AV-  and  V^ 
show  excellent  agreement  with  observed  erosion  behavior  of  the 
specimen. 
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A  'Similar  analysis  for  test  PY-ll  (Fig.  1?) 
indicated  that  the  mass  of  material  corresponding  to  AV^  and 
is  218  rag.  and  328  mg.,  respectively.  Microscopic  examina¬ 
tion  of  specimen  PY-ll  revealed  that  relatively  few  collisions 
were  occurring  over  th-  inboard  half  of  the  specimen.  When 
the  test  was  terminated,  the  surface  layer  at  the  outboard 
end  of  the  specimen  was  almost  completely  removed,  but  about 
one-sixth  of  the  specimen  surface  was  sparsely  eroded  at  the 
inboard  end.  A  moderately  dense  array  of  frustrums  covered 
the  central  portion  of  the  specimen.  Since  several  erosion 
mechanisms  were  operative  across  the  lateral  dimensions  of  the 
specimen,  it  is  diffi<  ult  to  correlate  the  anticipated  ero¬ 
sion  behavior  with  the  form  of  the  curve  in  Fig.  17,  However, 
on  the  basis  of  the  condition  of  the  eroded  specimen  surface 
at  the  end' of  the  test,  the  calculated  ’"asses  for  the  surface 
layer  are  quite  reasonable. 


The  weight- loss  data  for  the  70-raicron  glass 
beads  impacting  glass  specimens  PY-22  at  375  fpo  and  PY-23  at 
450  fps  and  a  bead  flew  rate  of  950  g/min  are  given  in  Fig.  18 
and  19*  A  second  glass  specimen  was  run  in  each  case  (PY-25 
and  PY~ 26)  to  obtain  the  weight  loss  during  the  «ar.ly  stages 
of  the  erosion  process  and  thesr  results  are  also  shown  in 
the  figures.  The  duration  of  the  initial  stage  of  the  erosion 
process  is  seen  to  vary  considerably  between  the  specimens 
tested*  however,  the  indication  in  Fig.  19  is  that  the  erosion 
rates  are  comparable  in  the  later  stages  of  the  erosion  pro¬ 
cess,  even  thougu  the  magnitude  of  the  weight  loss  may  not 
be  close  for  the  same  exposure  time. 


Again,  using  the  data  in  Table  VIII,  the 
mass  of  the  surface  layer  for  specimen  PY-22  and  PY-23  is 
11,9  mg  and  14V3  mg,  respectively.  The  surface  area  exposed 
to  particle  impacts  is  0.821  sq.ln.  Referring  to  Fig.  18  and 
19»  it  is  seen  that  the  70  micron  beads  do  not  produce  the 
bulge  in  the  weight-loss  data  as  in  the  case  of  the  290  micron 
beads,  The  steady- state  erosion  rate  is  quite  uniform  at 
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Weight- loss  data  for  speeimer  *:?- II  .impacted  hy 
2$0-micror<  glass  beads  at  37?  fp® 


an  impact,  velocity  of  &50  fps  but  shows  less  uniformity  at 
375  fps.  After  the  initial  period  depicted  in  Fig,  4  and  5 
the  process  of  material  removal  is  oy  fine-scale  chipping  or 
the  bulk  material. 

The  conical  f rust  rums  formed  on  the  surface 
oT  the  glass  specimens  by  70-micron  beads  impacting  at  375 
fps  are  shown  in  Fig.  20.  The  removal  of  material  is  in  the 
immediate  vicinity  of  each  ring  shortly  after  the  rings 
appear  on  the  surface,  Fig*  20  can  be  compared  with  Fig.  6 
where  a  dense  array  of  ring- fracture  •'  develops  before  large- 
scale  material  removal  takes  place.  The  condition  of  the 
surface  in  the  latter  stages  of  the  erosion  process  is  illus¬ 
trated  in  Fig.  21.  It  shows  the  uniform  fine-scale  chipping 
of  the  highly fractured. :  surface. .  The  scale  of  the  erosion 
particles  can  be  compared  with  those  due  to.  impacts  by  the 
290-micron  beads  as  shown  in  Fig.  6. 

.  3.3:  Oblique  impact  by  .Clase .Beads 

In  order  to  evaluate  the  extent  of  the  e rosier, 
damage,  produced  by  the  70-  and  290-mi. 'em  glass  beads  and- the. 
90-micron  sand  impact  inf  obliquely,  as  compared  direct 
impacts*  the  'partiolevfiow  ratej  -and -the  particle  vel-eciiv 
normal  to  the  specimen's  d.ur^aeKhere.  main*-*’  *»  v ne  values 
.used  for  .-direct  impacts .  thvde.r  tb»  total 

’ number  of  partible . Impacts by  the  .ap»ci.m?hs 
each  case  will  be  the  ^n^-s  41; ter  .equal  times  it 

hj  (oblique)  x'  (RFMV,  •■*  r  {-.-Iree*  lx  <9**\ 

where  $T  ir- -  f  pu**'  1  •-  t’.le  '■'specir-in  -experl- 

e'r^'.-..-- rr’ r^Vviur: .  ;v  U'  rotating  aim.  introducing-  $q*  (1 ) 

is'-v-  ?q. 


Figure  21  Eroded  surface  of  a  glass  specimen  at 

an  advanced,  stage  of  the  erosion  process 
due  to  70-micron  beads  impacting  at 
375  fps  (1167). 
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where  the  slight  difference  in  the  nozzle  dimension,  i,  for 
each  case  has  been  neglected  and  n  is  the  number  of  revolutions 
per  minute  for  direct  impacts.  It  is  seen  from  Eq.  (26)  that 
the  weight-loss  data  for  oblique  impacts  can  be  compared  with 
that  for  direct  impacts  by  multiplying  the  values  of  the 
weight  loss  for  the  oblique  impacts  by  the  ratio  of  the  exposed 
areas.  For  the  290  micron  beads  the  exposed  area  for  direct 
impacts  is  2.25  sq.in.  while  it  is  0.82  sq.in,  for  the  TO 
micron  beads  using  the  experimental  arrangement  shown  in 
Fig.  1.  The  exposed  area  for  oblique  impacts  is  .1.81  sq.  in. 
in  both  cases.  It  follows  from  Eq. (26)  that 
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The  ratio  of  the  areas  for  the '2’90  micron  beads  is  1.25  and 
0.453  for  the  70  micron  beads.  The  weight-loss  data  for 
oblique  impacts  are  plotted  both  In  terms  of  the  actual  mea¬ 
surements  and  with  the  adjustment  for  the  difference  in  the 
eroding  areas. 


Oblique  Impact  by  290-Micron  Glass  Beads 

The  weight- loss  data  for  normal  impact  velo¬ 
cities  of  200  and  300  fps  and  a  bead  flow  rate  of  370  g/min. 
are  given  in  Fig.  22  and  23.  It  is  noted  that  the  convex 
hump  appearing  in  the  weight-loss  data  for  direct  impacts 
(Fig.  16  and  17)  after  the  initiation  stage  corresponding  to 
the  sequence  of  erosion  mechanisms  described  in  Section  3.2.1 
is  also  present  for  the  oblique  impacts.  Comparison  of  the 
curves  for  a  normal  impact  velocity  of  200  fps  shows  that  the 
weight  loss  for  the  oblique  impacts  is  from  two  to  three  times 
higher  than  that  for  direct  impacts  at  the  same  values  of  the 
exposure  time.  The  magnitude  of  the  weight  loss  can  be  brought 


into  coincidence  with  that  In  Pig.  16  if  the  time  scale  in 
Fig.  22  is  expanded  by  a  factor  of  two.,. 

Microscopic  examination  of  a  series  of  incre¬ 
mentally  eroded  specimens  revealed  that  the  operative  erosion 
mechanisms  were  identical  to  those  already- discovered  for 
direct  impacts  of  290  micron  beads.  At  the  termination  of 
test  PY-J+  material  was  completely  removed  from  the  inter¬ 
stitial  regions  and  chipping  of  the  protruding  frustrums  was 
taking  place.  Test  PY-35  was  terminated  at  a  point  where  the 
exposed  frustrums  were  still  intact  and  a  small  portion  of 
interstitial  material  could  still  be  observed  on  the  eroded 
surface.  These  observations  correspond  to  the  characteristic 
shapes  of  the  weight-loss  data  reported  in  Fig.  22  and  23. 

The  particular  phase  of  the  overall  erosion  process  pred-omin- 
ating  as  the  exposure  time  increased  appeared  uniformly  over 
the  face  of  the  specimen. 

Detailed  microscopic  examination  of  the  cone 
fracture  showed  that  they  were  circular  at  the  surface  of  the 
specimen.  Histograms  of  the  diameter  of  the  surface  ring 
fracture  are  given  in  Fig.  24  for  specimen  PY-38  and  PY-39 
and  in  Fig. ■ 25  for  specimen  PY-36.  The  mean  diameter  of  the 
ring  fracture  is  3.28  mils  at  200  fps  and  3*80  mils  at  375  fps. 
These  values  are  based  on  a  sampling  of  more  than  125  ring 
fractures.  The  mean  cone  penetration  depth  below  the  surface 
of  the  specimen  is  1  mil  (25ft)'- at  200  fps  and  1.6  mil  (40p) 
at  375  fps.  These  depths  are  somewhat  below  those  computed 
in  Table  VIII.  The  corresponding  cone  fracture  data  for  direct 
impacts  are  a  mean  ring  diameter  of  3.17  mils  and  cone  depth 
of  25- 30m  at  200  fps  and  a  mean  ring  diameter  of  3*72  mils 
and  cone  depth  of  65g  at  375  fps. 


fracture  diameters  on  the  surface 


3.3.2  Oblique  Impact  by  70-Micron  Glass  Beads 

Weight-loss  data  were  obtained  at  normal 
impact  velocities  of  300  and  375  fps  at  the  same  bead  flow 
rate  of  950  g/min.  used  in  the  direct  impact  tests.  Fig.  26 
and  27  indicate  that  a  steady-state  erosion  rate  has  not.  been 
achieved  even  at  exposure  times  which  are  well  in  excess  of 
those  required  to  remove  the  surface  layer  of  material  from 
the  specimen.  Comparison  of  the  weight-loss  data  obtained 
for  oblique  and  direct  impacts  at  a  normal  velocity  of  375 
fps  (Fig.  18  and  27)  indicates  that  the  values  of  the  weight 
loss  for  glass  plate  at  an  impact  angle  of  50°  are  from  two 
to  four  times  greater  than  the  direct ' impact  case.  Again, 
as  in  the  case  of  the  290-micron  beads,  closer  agreement  can 
be  achieved  if  the  time  axis  in  Fig.  27  is  expanded  by  a  fac¬ 
tor  of  two.  From  the  limited,  data  available  it  appears  that 
the  erosion  rates  are  comparable  at  the  same  exposure  times, 
however  the  probability  for  ring  fracture  formation  is  greatly 
enhanced  for  the  oblique  impacts. 

.  •  .  ■  •  ■  -  •  '  •  •  .  v 

3.4  Direct  and:  Oblique  Impact  of  Irregularly-Shaped 
Particles 

Most  of  the  work  on  solid  particle  erosion  has  been 
involved .with  the  long-time,  steady-state  erosion  rates. 

Models  describing  the  steady-state  erosion  by  solid  particles 
have  been  developed  on  the  basis  of  laboratory  erosion  tests 
employing  sandblast  testers.  The  basic  premise  employed  in 
these  Semi-empirical  models  is  that  a  certain  amount  of 
energy  is  required  to  remove  material  from  the  surface  of  a 
target  material  and  that  each  material  has  a  characteristic 
erosion  resistance.  Two  different  modes  of  particulate  ero¬ 
sion  are  distinguished  empirically  corresponding  to  ductile 
and  brittle  erosion  behavior  of  the  target  material.  In 
sandblast,  experiments  uniformly-graded  sand  particles  carried 
by  a  high-velocity  airst  re  am  Impact,  the  target  material.  The 
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result's  of  the  test  are  reported  in  terms  of  the  mass  of 
material  removed  per  mass  of  impacting  particles. 

The  research  reported  here  is  an  initial  effort  in 
extending  the  concepts  and  experimental  correlations  esta¬ 
blished  for  spherical  particles  to  irregularly-shaped  quartz 
sand  particles.  The  experimental  procedure  employing  the  AFML- 
Bell  rotating  arm  differs  from  the  solid  particle  erosion 
tests  conducted  in  the  sandblast  apparatus  in  that  the  com¬ 
ponent  of  the  impact  velocity  normal  to  the  specimen  surface 
was  held  constant  in  the  direct  and  oblique  impact  tests. 

These  conditions  afforded  an  opportunity  to  evaluate  the 
influence  of  the  velocity  component  tangent  to  the  surface  of 
the  specimen  on  the  erosion  behavior  of  the  glass  plate. 

Orly  a  slight  increase  in  erosion  rate  was  found  and  no  real 
distinction  could  be  observed  microscopically  in  the  damage 
due  to  obliquely  impacting  irregularly- shaped  particles 
compared  with  direct  impacts. 

•  . 

3.^.1  Erosion  of  Materials  by  Irregularly- Shaped 

Solid  Particles 

Finnie^’ developed  a  mathematical  model 
for  describing  the  erosion  of  ductile  materials  as  a  func¬ 
tion  of  impingement  angle  based  on  an  erosion  mechanism  whereby 
material  removal  is  entirely  due  to  the  displacing  or  cutting 
action  of  an  abrasive  particle  impacting  a  perfectly  plastic 
material  surface.  The  equation  derived  by  Finnie  was  moder¬ 
ately  successful  in  describing  the  angular  dependence  of  the 
steady-state  erosion  rate  for  a  given  set  of  particle  and 
regime  parameters  on  the  basis  of  a  single  erosion  test.  The 
model  predicted  that  there  would  be  no  erosion  of  ductile 
metals  at  normal  impingement  angles,  however  this  result  does 
not  agree  with  the  experimental  findings.  More  general  models 
were  proposed  by  Bitter  and  Neilson  and  Gilchrist ^ ^ 
which  included  both  cutting  wear  (proposed  by  Finnie)  and 
erosion  due  to  repeated  deformation.  Repeated  deformation  is 
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responsible  for  cracking  the  surface  of  a  brittle  material 
with  subsequent  intersection  of  cracks  and  material  removal . 
These  equations  provide  better  curve  fit' tine  of  the  experi¬ 
mental  data.  These  formulas  include  the  idea  of  a  particle 
energy  threshold  below  which  deformation  erosion  ceases  and 
a  minimum  effective  angle  of  impingement  below  which  ductile 
erosion  ceases.  Implementation  of  Bitter's  model  requires 
knowledge  of  the  energies  required  to  remove  a  unit  volume 
of  a  material  by  cutting  wear  and  deformation  wear,  respec¬ 
tively.  parameters  are  determined  by  performing  at 

least  two  tests  and  back-calculating  to  determine  the  neces¬ 
sary  values.  Neilson  and  Gilchrist^1'  provided  simpler 
equations  for  deformation  wear  than  those  of  Bitter^'  which 
were  able  to  account  for  the  angular  dependence  of  the  erosion 
process  and  included  another  parameter  to  be  evaluated 
mentally  which  is  the  normal  velocity  component  below  which 
no  erosion  takes  place.  ,  These  models  have  no  predictive  capa¬ 
city  to  describe  the  steady- state  erosi on  of  brittle 
ductile  materials  eroded  by  natural  contaminants  caused  ■-•?? f 
a  distribution  of  particle  sites  and  shapes. 

In  general*,  the  erosion  rate  is  a  function, 
of  a  number  of  parameters  associated  with  the  Impact,  regime, 
the  impacting  particles,  the  target,  and  the  part tele/t arget 
interaction.  The  accumulation  of  experimental  data,  has  shewn 
the  relative  importance  of  acme  of  the  impact  regime  para¬ 
meters  which  inc iude  the  distribut ion  of  t he  part isle  impact 
veloeiti  impingement  angle,  and  concentration.  Tt*  par¬ 
ticle  parameters  are  the  distribution  of  particle  aikes, 
shapes,  and  quantitative  measures  of  physical  and  mechanics) 
properties*  The  target,  parameters  include  quart  it  a? ive  mea¬ 
sures  of  its  physical  end  mechanical  proper4 ies  and  the 
condition  of  the  surface  before  exposure  to  the  eroding  envi¬ 
ronment. 


The  experimental  data  obtained  for  glass 
bead  impacts  on  glass  places  have  successfully  defined  the 
operative  erosion  mechanisms  associated  with  a  set  of  con¬ 
ditions  which  could  be  used  to  gain  understanding  for  model¬ 
ing  the  general  erosion  process  for  arbitrarily- shaped  quartz 
sand  particles.  The  use  of  glass  beads  and  glass  plates  made 
possible  the  introduction  of  mathematical  analysis  to  provide 
a  rational  basis  for  the  inclusion  of  particular  particle, 
target,  and  regime  parameters.  This  includes  the  influence 
of  particle  size  and  velocity;  however  further  extension  of 
the  model  to  sand  erosion  of  brittle  materials  requires 
detailed  information  on  the  quantitative  characterization  of 
irregularly- shaped  particles,  the  contribution  of  a  particular 
size  particle  in  a  distribution  of  particle  sizes  to  the  over¬ 
all  erosion  rate,  and  tne ' inf luence  of  the  angle  of  attack  on 
erosion  behavior. 


The  approximate  particle  size  distributions 
used  to  establish  some  correspondence  between  the  quartz  sand 
and  spherical  glass  beads  used  in  the  controlled  erosion  tests 
is  not  adequate  to  provide  a  quantitative  measure  of  the  geo¬ 
metry  of  irregularly- shaped  solid  particles.  In  order  to 
evaluate  the  relative  erosiveness  of  a  given  solid  particle 
environment  it  is  important  to  be  able  to  characterize  a  given 
sample  of  particulate  matter  by  means  of  simple  parameters 
which  define  the  shape  of  the  individual  particles.  A  survey 
of  the  literature  on  classification  of  particulate  matter 
which  is  of  general  interest  in  erosion  research  revealed  that 
there  is  no  adequate  set  of  parameters  to  provide  a  quantita¬ 
tive  description  of  particle  shape. 

(IP ) 

Hausnerv  *'  has  enumerated  several  require¬ 
ments  for  a  suitable  method  for  describing  the  shape  of  a 
particle.  These  requirements  and  assumptions  are  listed 
below  as  amended  by  Heywood^^  and  further  modified  here. 
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(1)  The  quantification  of  particle  shape 

is  a  three-dimensional  problem,  however 
from  the  standpoint  of  practicality  it 
is  treated  as  a  two-dimensional  problem. 

(2)  The  diameter,  or  any  other  quantity 
related  to  the  actual  size  of  the  par¬ 
ticle,  should  not  be  used  in  the 
characterization  of  particle  shape. 

,  Shape  >  is  a  geometrical  property  distinct 
''  from  size. 

(3)  "A  characterisation  of  -  particle  shape  •  ' 
;  must  consist  cf  a  series  of  factor's, 

the  knowledge  of  which  should  permit 
drawing  an  approximate  sketch  of  the 
particle. 

An  additional  observation  by  Hausner  is  that  even  the  best  . 
characterization  of  particle  shape  will  be  an  approximation. 

He  then  proceeds  to  define  three  parameters  which  are  descrip¬ 
tors  of  particle  shape:  the  elongation,  bulkiness,  and  surface 
factors.  These  factors  are  not  independent  of  those  proposed 
by  other  investigations. 

The  lack  of  uniqueness  in  specifying  the  size 
of  an  irregularly-shaped  particle  is  discussed  in  detail  by 
Kaye  especially  for  particle  with  re-entrant  particle  pro¬ 
files.  Statistically-evaluated  Feret's  and  Martin's  diameters 
are  compared  for  this  case.  Feret's  diameter  is  the  projected 
length  of  the  particle  profile  with  respect  to  a  fixed  direc¬ 
tion.  Martin's  diameter  is  the  length  of  the  line  through 
the  particle  profile  parallel  to  a  fixed  direction  which 
divides  the  particle  profile  into  two  equal  areas.  He  states 
that  Feret's  diameter,  which  is  consistent  with  the  shape  fac¬ 
tors  proposed  by  Hausner  and  others,  is  identified  as  the 
diameter  of  a  circle  with  a  perimeter  equivalent  to  that  of 
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the  particles,  whereas  Martin's  diameter  does  not  appear  to 
have  a  simple  relationship  to  any  obvious  particle  physical 
parameter. 


A  small  portion  of  the  quartz  sand  particles 
shown  in  Fig.  3  exhibit  re-entrant  particle  profiles.  Kaye^^ 
defines  a  convex  hull  which  does  not  include  convex  portions 
of  the  particle  profile  and  shows  how  Feret's  diameter  is 
still  meaningful  based  on  the  concept  of  a  convex  hull.  The 
involution,  of  the  particle  profile  would  be  the  difference 
between  the  actual  particle  perimeter  and  the  perimeter  of 
the  convex  hull.  In  general,  the  majority  of  the  quartz  sand 
particles  used  in  this  study  are  angular  with  a  small  per¬ 
centage  of  aeicular  particles:  the  percentage  of  acicular 
particles  increases  as  the  particle  size  decreases. 

The  particle  characteristics  usually  reported 
in  the  literature  on  solid  particle  erosion  are  the  mesh-size 
of  the  particles  used.  Only  Head  and  coworkers  (-^-l? )  (lave 
included  more  details  pertaining  to  particle  characteristics 
in  their  erosion  tests  and  modeling  studies. 

After  completion  of  an  extensive  testing 
program^1^,  Head  and  Harr^1^  attempted  to  develop  a  model 
applicable  to  natural  soils  which  includes  consideration  of 
the  energy  transferred  from  the  impinging  particles  to  the 
target,  the  natm  e  of  the  response  of  the  target,  and  the 
nature  of  the  erosive  agent  including  pertinent  descriptors 
of  composition,  angularity,  hardness,  and  size  distribution. 

The  result  of  this  effort  was  the  development  of  a  weighted 
regression  analysis  program  (WRAP  model)  in  which 

2 

A  =  f  Gg  ,  R,  sina,  cosa,  (28) 

where  A  is  the  erosion  rate  measured  in  volume  removed  per 
unit  mass  of  dust,  V  is  the  effective  velocity  of  the  mix,  R 


is  the  effective  roundness  derived  from  a  quantitative  evalu¬ 
ation  of  the  particle  shapes,  H  is  the  effective  hardness  of 
the  dust  particles,  a  is  the  particle  impingement  angle,  B 
is  the  hardness  of  the  target,  and  E  is  the  erosion  resistance 
per  unit  volume  of  target.  Head  and  Harr  use  the  effective 
parameters  which  they  define  as  the  sum  of  the  parameters 
associated  with  the  individual  size  components  of  the  dust 
weighted  on  the  basis  of  the  grain-size  distribution  of  the 
dust.  A  different  form. of  interaction  between  the  dependent 
variables  was  found  for  brittle  and  ductile  materials. 

The  WRAP  model  includes  the  influence  of 
particle  size  as  a  weighting  function  for  other  parameters 
and  a  single  factor,  Wadell*. s  roundness  parameter^® \  is 
employed  to  characterize  the  shape  of  the  particles.  On  the 
basis  of  the  earlier  discussion,  this  mode  of  particles  char¬ 
acterization  is  probably  not  adequate  to  provide  good  pre¬ 
dictions  of  the  erosiveness  of  a  range  of  particle  geometries 
and  size  distributions.  An  alternative  statistical  formula¬ 
tion  will  be  considered. 

The  predictive  ability  of  the  WRAP  model  for 
brittle  materials  was  not  determined  by  Head  and  Harr^^, 
however  the  predictive  capability  of  the  WRAP  model  was  sub¬ 
sequently  evaluated  for  two  ductile  metals  (type  302  stain¬ 
less  steel  and  6061-T6  aluminum)  eroded  by  220  mesh  alumina 
particles  and  220  mesh  fluorite  particles ^ \  When  the 
model  was  adjusted  for  a  higher  maximum  particle  velocity, 

300  m/sec.  as  opposed  to  185  m/sec  in  the  earlier  work  imparted 
by  the  carrier  gas  and  changes  in  the  dimensions  of  the  blast 
tube,  agreement  between  the  model  predictions  and  the  experi¬ 
mental  data  for  alumina  particles  is  claimed  for  impingement 
angles  between  20°  and  60°  with  good  estimates  of  the  maxi¬ 
mum  erosion  rate;  however,  the  general  trend  of  the  model 
prediction  and  the  experimental  data  as  a  function  of  impinge¬ 
ment.  angle  show  little  resemblance  to  each  other.  The 
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comparison  is  less  encouraging  for  the  fluorite  particles 
where  the  pred i cl  ed  maximum  erosion  -■■ates  underestimate  the 
actual  values  by  a  fact  or  of  two  for  the  stainless  steel  and 
a  factor  of  four  for  the  aluminum  alloy.  It  is  interesting 
to  note  that  the  fluori'-'  particles  with  a  Moh  hardness  of 
4.0  and  higher  particle  roundness  than  the  alumina  were  sig¬ 
nificantly  more  erosive  than  tlv  alumina  particles  with  a 
Moh  hardness  of  9.0.  In  the  case  of  f>06l-T6  aluminum  there 
is  nearly  an  order  of  magnitude  difference  in  the  measured 
erosion  rates  for  these  two  particles.  The  WRAP  model  predict 
that  the  erosion  rate  increases  with  pari icle  hardness  and 
roundness.  Evidently,  there  are  other  factors  which  must  be 
taken  into  account  with  regard  to  particle  characterization 
and  testing  procedures. 

3.4.2  Weight- loss  Data 

The  weight- loss  data  for  oblique  impacts 
were  acquired  under  conditions  which  differ  from  those 
reported  for  sandblast  experiments.  Tn  the  sandblast  tests 
the  velocity  of  th°  airstream  is  held  constant  and  the  inclin¬ 
ation  of  the  specimen  with  the  airstream  is  changed.  Under 
these  conditions  the  velocity  component  normal  to  the  surface 
of  the  plate  varies  as  the  cosine  of  the  angle  between  the 
normal  to  the  plane  face  of  the  specimen  and  the  direction 
of  the  air-flow.  The  experiment  s  conducted  in  the  AFML-Bell 
rotating  arm  facility  were  designed  so  that  the  particle 
Impact  velocity  normal  to  the  plane  face  of  the  glass  speci¬ 
men  was  held  constant.  This  approach  would  allow  a  better- 
evaluation  of  the  change  in  the  effective  erosiveness  of  the 
particles  with  angle  of  attack  without  introducing  «!dltioral 
nonlinear  effects  associated  with  a  variable  Impact  velocity. 
Due  to  experimental  difficulties  encountered  in  obtaining 
data  in  the  runs  for  glass  bead  impacts,  a  comparable  set  of 
data  is  not  available  for  both  direct  and  oblique  impacts. 

It  is  found  on  the  basis  of  the  runs  which  can  be  compared 


that  the  steady-state  erosion  rates  are  the  same  to  within 
experimental  error  for  the  glass  beads  impacting  at  an  angle 
of  50°  and  normal  to  the  plane  surface  of  the  specimen.  Data 
are  more  complete  for  the  90-micron  quartz  sand.  To  within 
the  degree  of  variability  found  for  the  erosion  behavior  of 
glass  specimens,  the  steady-state  erosion  rates  can  be  con¬ 
sidered  equivalent  with  the  oblique  impacts  tending  to  be  just 
slightly  more  erosive  than  the  direct  impacts. 


3. 4.2.1  Direct  Impact  by  90-Micron  Sand 

From  the  sieve  analysis  reported  in  Fig.  2 
the  average  lateral  dimension  of  the  quartz  sand  was  found  to 
be  roughly  90  microns.  The  flow  rate  for  the  sand  particles 
had  to  be  reduced  to  260  g/min.  from  the  950  g/min.  used  with 
the  glass  beads,  since  the  weight  loss  for  sand  was  greatly 
accelerated  over  that  observed  for  the  70-micron  •»Tass  beads. 

The  weight- loss  data  at  impact  velocities  of  200,  300,  and 
375  fps  are  provided  in  Fig.  28  to  30. 

In  comparing  the  weight-loss  data  for  the 
sand  particles  with  that  obtained  using  glass  beads,  it  is 
necessary  to  take  the  different  bead  flow  rates  into  account. 

In  the  present  context  the  number  of  impacts  per  revolution 
recorded  in  Table  V  would  provide  a  better  basis  for  compari¬ 
son.  Even  after  accounting  for  these  differences  the  experi¬ 
mental  conditions,  the  sand  particles  are  found  to  be  consider¬ 
ably  more  erosive  than  the  glass  beads. 

The  general  character  of  the  weight- loss 
curves  shows  that  after  an  initial  stage  of  unsteady  erosion 
rate,  due  in  part  to  the  fact  that  the  weight  measurements 
are  only  accurate  to  1  mg  and  a  probability  f.-r  measuring 
errors  exists,  a  uniform  rate  of  material  removal  Is  established 
at  each  velocity.  The  values  for  the  uniform  erosion  rates 
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fps.  (Test  PY-ifO). 


PF  30'  less,  for  direct  impacts  of  90-micrrn  sand  on  glass 

at  i?5  fps.  (Test  rY-U2) 


are  plotted  as  a  function  of  velocity  . in  Fig.  31.  The  raw 
data  can  be  described  by  the  following  relation, 

IE  =  (3.35  x  10’ 10 )V3* 65  (29) 

where  I  is  the  erosion  rate  expressed  in  mg/sec  and  V  is  the 

fcj 

impact  velocity  in  fps. 

i 

3.^.  2.2  Oblique  Impact  by  90-Micron  Sand 

The  sand  flow  rate. was  maintained  at  the 
same  rate  that  was  us^d  for  the  direct  impacts  which  was 
260  g/min.  Microscopic  examination  of  the  eroded  surfaces  did 
not  reveal  any  differences  between  the  damage  due  to  oblique 
impacts  and  that  associated  with  direct  impacts. 

The  weight- loss  data  at  normal  impact 
velocities  of  200,  300  and  3?5  fpa  are  shown  in  Fig.  32  to 
.  3Uv  After  a. relatively  brief  initial  period  of  unsteady  mate¬ 
rial  removal  a  steady  rate  of  material  removal  Is  established. 
The  steady-state  erosicrv  rates  at  each  velocity  are  recorded 
in  Table  IX.  The  data  at  3OO  fps  appears  to  be  somewhat 
anomalous,  however  the  data  at  the  other  two  velocities  are 
consistent  Indicating  that  the  erosion  rates  are  slightly 
higher  for  oblique  impacts  at  an  angle  of  attack  of  HO''  to 
the  surface  of  the  specimens  than  direct  impacts  when  the 
normal  impact  velocity  is  held  constant. 

3.^.3  Microscopic  Examination 

Microscopic  examination  of  specimens  FY-U3 
*0  U6  and  to  55  for  the  direct  and  oblique  impact  of 

the  90-micron  sand  particles  under  the  conditions  listed  i«, 
Table  I  indicated  that  there  was  very  little  difference  in 
the  form  of  the  erosion  damage  in  each  case.  A  series  of 


Weight-Lot*  Rate  (mg/iecJ 


TABLE  IX 

STEADY  RATE  OF  MATERIAL  REMOVAL  FROM  GLASS 
SPECIMENS  FOR  DIRECT  AND  OBLIQUE  IMPACT  OF 
90-MICRON  SAND  AT  VARIOUS  PARTICLE  IMPACT  VELOCITIES 


Normal 


Iropac  t 
Velocity 
(fps) 


Erosion  Rale  (nig/ce0 /unit  area) 


Impac  t 
Angle 
0° 


Impact 

Angle 

50° 


200 


0.100 


0.105 


300  0.425  0.224 

375  1.00  1.26 
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micrographs  is  provided  in  Fig.  35  and  36  which  shows  repre¬ 
sentative  damage  areas  on  the  surface  of  c aeh  specimen  and 
comparisons  between  equivalent  impact  conditions  for  direct 
and  oblique  impacts.  These  areas  were  selec-ed  in  a  random 
manner,  however  a  qualitative  evaluation  of  the  resulting 
damage  in  each  instance  shows  minor  differences  in  the  impact 
indentations  for  direct  and  oblique  impacts.  The  impacts  at 
375  fps  are  the  most  damaging,  and  this  effect  appears  in  the 
larger  indentations  which  are  formed  at  this  velocity.  Quan¬ 
titative  evaluation  of  the  damage  is  required  Ln  the  modeling 
studies  in  terms  of  specific  characteristics  of  the  particle 
indentation:  the  distribution  of  size,  shape  and  penetration 
as  a  function  of  impact  velocity  and  exposure  time.  This 
effort  parallels  that  of  characterizing  the  sand  particles 
comprising  the  erosive  medium  and  will  be  investigated  in  the 
near  future. 


The  progression  of  damage  on  the  eroded  sur¬ 
face  of  specimen  PY-50  can  be  seen  over  a  limited  range  of 
the  exposure  time  in  Fig.  37.  The  same  area  is  monitored  for 
each  time  increment.  Once  a  localized  pit  develops  it  does 
tend  to  grow;  however  the  nucleation  rate  is  high  for  the 
sand  particles,  so  numerous  additional  pits  v. ill  be  formed 
in  the  vicinity  of  the  first  one  before  it  has  grown  appreci¬ 
ably,  However,  there  is  still  abundant  evidence  that  pit 
growth  is  preferential. 


3.5  Liquid  Droplet  Impacts 

A  number  of  mathematical  analyses  applicable  to  the 
quantitative  evaluation  of  the  transient  stress  conditions 
which  prevail  within  a  material  body  when  it  is  impacted  by 
a  spherical  liquid  drop  ere  considered  in  this  section. 

While  the  complete  mathematical  problem  for  a  drop  striking 
a  deformable  material  surface  at  moderate  velocities  is  beyond 
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the  existing  state-of-the-art,  a  number  of  reasonable 
approximations  to  the  general  problem  are  presented  here. 

Only  a  limited  amount  of  numerical  results  are  provided  due 
to  the  complex  nature  of  the  mathematical  solutions  which  are 
available.  This  research  is  only  an  interim  report  on  the 
advances  in  computer  modeling  which  offer  considerable  poten¬ 
tial  for  evaluation  and  development  of  material  systems  with 
greatly  enhanced  erosion  resistance.  The  required  computer 
codes  will  be  compiled  and  developed  further  to  achieve  this 
objective  in  future  work  on  this  program. 

The  erosion  damage  mechanism  has  to  be  established 
for  different  types  of  material  response.  For  brittle  mate¬ 
rials  the  magnitude  and  duration  of  the  critical  radial  sur¬ 
face  stress  component  may  be  the  most  important  factor.  In  a 
ductile  material  the  critical  shear  stress  trajectories 
would  be  more  significant.  The  availability  of  data  on  tran¬ 
sient  stress  distributions  within  the  interior  of  the  solid 
body  provides  a  means  for  quantitative  evaluation  of  the 
critical  conditions  under  which  failure  of  a  particular  mate¬ 
ria."  system  is  likely  to  occur  and  to  evaluate  the  erosion 
resistance  potential  for  materials  from  computer  modeling 
studies. 

Table  I  lists  the  tests  conducted  on  glass  specimens 
exposed  to  a  standard  rainfall  of  1  in/hr.  of  1.8  mm  drops 
on  the  AFML-Bell  rotating  arm.  Although  the  glass  specimens 
impacted  by  glass  beads  could  sustain  velocities  around  500 
fps  without  fracturing,  this  was  not  the  case  in  the  rain 
erosion  tests.  It  was  found  that  rain  erosion  damage  on 
glass  was  inconsequential  below  400  fps.  In  order  to  relieve 
the  stress  applied  to  the  specimen  ir,  the  specimen  holder  at. 
the  end  of  the  rotating  arm  at  velocities  above  400  fps,  the 
glass  specimen  was  shortened  and  the  specimen  then  bonded  to 
an  aluminum  backing  plate  to  reduce  problems  with  fracturing 
of  the  ends  of  the  glass  specimens.  This  arrangement  was 


adequate  to  act. 'eve  the  velocities  required  to  produce  erosion 
damage  by  water  particle  Impacts  on  the  glass  specimens  but 
not  for  the  exposure  times  required  for  the  acquisition  of 
weight-loss  data. 

The  number  of  impacts  that  the  specimen  receives  in 
the  rain  fields  produced  in  the  AFML-Bell  erosion  facility 
can  be  calculated  in  a  manner  similar  to  that  developed  for 
solid  particle  impacts.  The  r-ainfield  is  constant  o  rer  the 
complete  area  described  by  the  specimen  in  one  revolution. 

This  area  is 


A  =  Jr(r„2- r,2) 


(30) 


where  r^«108  in.,  r0*»lll  in.,  are  the  distances  from  the  cen¬ 
ter  of  rotation  to  the  inboard  and  outboard  ends  of  the  speci¬ 
men  on  the  rotating  arm.  The  number  of  drops  per  hour,  N,  is 
given  by  ;  . 

N  -  A  •  (51) 

?r  ' 

where  j  is  the  magnitude  of  the  ra infield  in  in, /hr.  and  r 
is  the  radius  oi  a  raindrop.  According  to  Eq. (22)  (Reference 
1)  the  density  of  a  given  flow  rate  of  particles  in  a  gravi¬ 
tational  field  at  a  distance  below  the  point  from  which  they 
were  released  with  zero  initial  velocity  is 


>(x) 


(3?) 


where  g*  32.2  ft. /sec.  and  X  is  measured  in  feet.  The 
plane  of  the  rotating  specimen  in  the  erosion  facility  is 
about  6  ft.  below  the  apex  of  the  trajectory  of  the  spray  of 
fluid  droplets.  To  obtain  the  number  of  droplets  the  speci¬ 
men  encounters  per  revolution,  the  density  from  Fq.  ( 3?)  Is 
multiplied  by  the  volume  the  specimen  passes  through  per 
revolution. 


99- 


Nd 


(33) 


43, 200>/5gx 


J(re2-ri2)^ 

57,600r^/5gx 


where  d  is  the  height  of  the  exposed  area  of  the  specimen's 
surface  in  Inches. 


For  the  aercdynamically- shaped  specimens  usually 
employed  in  erosion  tests  on  the  rotating  arm  and  exposed  to 
1  in. /hr.  rainfall  of  1.8  mm  drops,  d  =  0.5  in.,  so 

Nj  -  6.5  impact 3 /revolution 

The  flat-facer!  glass  specimens  exposed  to  the  same  rainfall 
encounter  9*?’:5  drops/r^volution.  These  estimates  of  the 
impacts  the  specimen  receives  will  be  used  in  evaluating  the 
erosion  behavior  of  the  materials  investigated  in  this  program. 


3*5. 1  Mechanics  of  Liquid  Prop  Collisions  With 
Deformable  Surfaces 

The  damage  produced  on  solid  surfaces  due  to 
liquid  particle  impingement  arises  from  two  sources:  the  high 
pressure  which  is  generated  over  the  expanding  impact  area  in 
order  to  bring  the  drop  to  a  sudden  stop  and  the  high-velocity 
lateral  flow  of  ttv  liquid  subsequently  escaping  from  t tv 
high-pressure  zone.  Detailed  knowledge  of  the  flow  patterns 
and  impact  pressures  developed  during  the  collision  of  a 
liquid  drop  and  a  deformable  plane  surface  is  still  unavail¬ 
able.  In  this  case  the  governing  field  equations  and  boundary 
conditions  are  coupled  and  pose  a  very  difficult  mathematical 
problem.  While  it  is  a  formidable  tack  to  solve  the  coupled 
mathematical  problem,  approximations  to  the  general  physical 
conditions  can  be  made.  Itv^se  alternative  approaches  fall 
into  three  major  categories; 


Evaluation  of  the  magnitude  of  the  pressure 
developed  at  the  interface  between  the  drop 
and  the  solid  surface  using  a  one-dimensional 
shock  wave  analysis  which  accounts  for  the 
compressibility  of  both  the  liquid  and  th» 
solid. 

Evaluation  of  the  pressures  and  deformations 
associated  with  the  drop  assuming  the  solid 
surface  is  rigid. 

Evaluation  of  the  stresses  and  deformations 
within  the  solid  after  prescribing  a  pres-, 
sure  distribution  which  is  assumed  repre- 
sent.-j.tt re  of  that  developed  at  the  interface 
between  the  drop  and  the  surface  of  .the 
solid  at  some  stage  of  the  collision  process. 

A  number  of  high-speed  photographic  studies'**’ 
provide  qualitative  insight  into  -he  phenomena  associated  with 
th?  single  impact  of  e  liquid  droplet  on  a  deformable  surface. 
Liquid  jets  striking  a  solid  surrace  at  ?hOO  fps  exhibited  a 
rise. 'time  of  IpS  for  the  peak  pressure  to  be  obtained  which 
iheri  decayed  to  a  much  lower  value  within  2  to  3ps  wtilch  is 
the  time  it  took  release  waves  to  propagate  into  the  interior 
©f  the  ,let.^0^  The  jet  continues  to  impact  the  surface  for 
another  ?Ops,  however  the  liquid/solid  interfacial  pressure  is 
at  a  greatly  reduced  level  assumed  to  be  equal  to  the  stagna¬ 
tion  pressure  along  the  axis  of  symmetry-'  The  experimental  data 
on  spherical  drops  striking  a  solid  surface  st  ill  cannot  pro¬ 
vide  a  quantitative  picture  of  the  temporal  magnitude  and  dis¬ 
tribution  of  the  liquid -’Solid  interfaclal  pressure. 

Detailed  analytical  studies  have  only  been 
provided  for  the  pressure  distribution  due  to  a  liquid  drop 
striking  nnooth,  rigid  surfaces. 5  f  Huang/*^  in  parti¬ 
cular  has  pres*  nfced  the  complete  mechanics  analysis  with 


computer  solutions  to  indicate  that  time-evolution  of  the 
impact  pressures  and  lateral  outflow  at  impact  velocities  of 
980  and  2^50  fps.  Huang's  calculations  support  Engel's  con¬ 
tention  .that  the  pressure  buildup  and  lateral  outflow 
occur  simultaneously  after  the  initial  contact  is  -made, 
although  the  jetting  is  not  appreciable  during  the  early  stage 
of  impact.  This  contradicts  the  more  widely  held  opinion  of 
Morris^8^,  Heymann^27^,  and  Bowden  and  Field ^  that'  lateral 
flow  cannot  begin  before  the  compression  wave  in  the  fluid 
moves  ahedd  of  the  circular  boundary  between  the  drop  and. the 
solid  surface.  At  approximately  30Q  fps,  these  calculations 

indicate  that  the  c  ri t  ical  edge  angle  is  on  the  order  of  a 

r  o-?  \ 

few  degrees.  Brunt  on  and  Camus’'-.,  however,  provide  photo-  •. 
graphic  evidence  that  under  these  conditions  the  critical 
angle  falls  in  the  range  of  from  13°  to  19°.*-  or  the  radius  of 
the  contact  area  at  this  point  .is  from  0. 2$  to  0,33  times. the 
radius  of  the  d**op»  These  results  arc  consistent  with  the 
t ime-dependeht  drop  corjf igurationa  computed  by  Huang  and  sup¬ 
ported  by  t?T§  .computer  analysis  of.  ■  >y 


Huang ‘ s  results ‘ *  *  ‘  also  indicate  thst 

the  ratio  of  the.  lateral  jetting  velocity  to  the  iropacl  Velo¬ 
city  decreases  as  the  impact  velocity  increases .  fhl.a  sum 
effect  is  evident  in  the  photographic  ^udy  of  .-.fenKlrs  and 
BooKtr>^'  where  to  3-1/2-min  drops  were  Impacted  „6n  .aynee't- 
su r faebfoyfr jg .  vsXoe t tyy  ‘f&fsga  of:  ’  <00  to.  fpc  V  -A  plot  of 
lateral  yclvsclty  as  ;a  function  of  thtj&apac*'  v?;loc.tt,y  shop* 

the  rati^  of  the  velocity  of  lateral pet :f loyz  ic  velo¬ 

city  of  impact  it  approximate!  6  at  an  impact  velocity'  of  1 
300  fps  while  it  is  only  2  at  250?  fps.*  fyall's •22-  mea¬ 
surement  of  radial  outflow  for  a  2  m  drop  striking  a  . plate 
Hf-  poiyfesthylacthae ryla  teat  approx itsateiV  1  0£*y  fps  K‘cu  1 1 
ih  roughly  the  same  ratio  between  the'  1  adial  joutflow  and  tne 
normal  impact  velocity  as  found  fey  -iters inland  Booker,  y 
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Brunt  on  and  Camus'  have  also  provided 

phonographic- -'e-vldifUleo  of  cavitation  bubbles  within  an 

Impact l;hg  vi,rc/pv  .cavitation  bubbles  *wr-re  found  at  the  liquid/ 

solid-  interface  -qurlniqthe  early  stage  of  impact  due  to  a 

negative  pressure-;  wtili'ch 'results  from,  the  expansion  of  the  com- 

pre  ssed;:  liquid  whe-tq  ^idbways  flow  begins.  Brunt  or.  and  Camus 

explain  that  the:- expatKUnCl  Iquiid  overshoots  and  a  negative 

pressure  fon?.s  at  ttie  surface  propagating  Into  the  liquid  as 

a  tension  tail  fg>Vbbe'  m&tn'  ct'mpressidn-.wave.  '.Cavitation 

buholeswefe  ^tecffoond  vbivm-jwar'  ♦■he  rear  face  of  the 

impacting  drop  due-'-hf-  'hhfr.BfeYX^e  Vlon.  of-  t  he-  impact  pressure 

wave  f rbro  .  t  fe’  ro$»-r  nurfp-*.:  of  J  he  (irco.  The  occurrence  of 

oavt-  a/  loif  bubble  in  Yfcis .vrstPi^rf-Wa^  originally  suggested,  by 
■  foiiV:  -V-  -  '  ' 

Ehfge.lv "  and  sUifoortbpi  hr  *  he  numerical  ealCUTa- 

f^j*  :  -  • 

t  iihis  of;  T:UAhffv'.  .(*n  -the' <SvYer  hantf  the  phriograchie.  . 
s f  ;..vd ! £$■ •  c c'is.dv it t  ed . by  1  ;  ..?/}*?■.  a  r  mm  .water  drop  Impact- 

/iffg  'a'.plaip.  ,;.%f  .^iymebv,;^e*haer;ylatV  af^OO  fps  did  hot 
e-thlbl  t  i irlis  e ff  ec t ;  Y  ffyai-i  '•'bS'pgt? . ’  fitted  t  hat  t he  pressur  & 
:%A«p.;feay'$*  afe;^\iat#cd.  ^S^t-n  vthe  '  dropl  ■. 


Y  .Thus,-  Huab:gVb  Yosspare  favorably 

tfiv-'b  t-ftsv.avY-1 1  sjbli  -phot-^f'  rapf  ic..  ;4  tu-i  lea  of  raind  rap  coll  is  ions » 
;lb  .gehe'rai^e  are-  inclined  vyo  Mil  eve  -that  Bumns's  analysis  •' 
is  a  ■featqtj&sl'y  Accurate,  at. Ion  of.  the.  dynamics  of  a 

ralr4rop;d*  ri>rbh:g  f,  ri^fd  |YrfaYe.v/  The  pressure  dlstr. Unit  ions 
whish  ee: 'proyfiej?  :4b  of- 1  and  spatial  coordinates 

b^ViSe?l:'in-  tpKanal>:^#a,  w«  wthtid  like  to  develop. 


.*  •••  :•/.  :  fV>ry  1:1  ttVe  r«nasvi*l  ty  can.  be  .found  ip  the 
'•I it 4iitikf*hln$ ;thg  nature  of  the  pressure  distribution  . 
arhich  .derM-bp*  _.<?■/&  solid  SMnface'  impacted  by  a  liquid  droplet, 
the  maghf  r  ide "  of  ifm'  pres  sure  put**  transmitted  to  a  defers;- 
Apia  m%ifrj*\.py  a  liquid  drop  baa  bfen  evaluated  in  temp 
cl  the  acoustic  impedance  of  t  be  liquid  suvd  tf$*  solid  material. 


0  C  V 
^000 


1  +  (  P0fco',plCl 1 


{3M 


where  and  are  the  respective  densities  of  the  liquid 

and  the  solid  and  C  and  C,  are  the  acoustic  velocities  in 

o  i 

the  liquid  and  solid,  respectively.  V  is  the  impact  velocity 
of  the  liquid  droplet.  This  approach  assumes  that  a  uniform"' 
pressure  is  impulsively  applied  over  the  entire  plane  surface 
of  a  semi-infinite  naif-space.  It  is  round  that  the  local 
impact  pressxires  in  raindrop  impacts  a‘  moderate  velocities 
are  of  a  magnitude  that  the  liquid  is  oomp re. ssed  . before  lateral 
outflow  is  initiated.  The  velocity  of  sound,  f\  will  be 
increased  due.  to  the  change  in  density .  Elementary  shock  wave 
theory  was  employed  Morris'  v ■  to.,  estimate  the  magnitude 
cf  the  initial  pressure  transmitted  to  s  sell;  under  these 
conditions., 

Huangs  numerical  apprcach  devised  .for 
obtaining  quant  Hat  lye  fo  n'.-iht'"  p-rvYsur#  distribution 

■associated  ■  wife,  a  conpr risible  .'liquid  dYop.  striding  a:.  rigid 
surface  indicates  that  ■ 

prv&mim* during,  the  \*my.  eari.y  ®tap*ayof  coilisHn. 

■'thr*  'may.tmvas  Value  of  the  is' reached  under  these  con¬ 

ditions.  •  *\s*  compression  of  the  drop  apntihpe&i.  the-dHlr  U>u-  ■" 
hi  on  beoer-es  nore-  unifom  vitfc'a  very  flight:  increase- 
in  .the  of  /the  pgsaaujr* , at  .thtf  ■.■periphery-'  of  ’he  con-., 

■tact'  area.'  iT;:.o  .ssa^hit udyy.ef.  the'.sssus.Sts.um  value  in-. the" pressure' 

dia?  riHU  Ion.  de^ryase:a:-.ai>'-  X9refai;  ssutnow-'becpmc^  morv 
d  bed.  nan  i,'.  . /-'ey  - .  ;.-v--  ■ 

/■'  .:- .  ■  .  .foniiderinn  ■ia^ui:yel.ntii'ios_  which'  are  less 

than H  !.-#  found,  analysis''  ?  '  prodHic  a 

.ra^iau^.  Pressure:  .c-f  . approx p „C  Jh  Y  'shl  le  Kuam?:''8 
iiu{reric'ai-.;'fe»u.l  t$V  .-'Indicate  -.shr  Y&xfmtus -pressure  decs  r.oi 
exceed.  0.  B  &  f-V  ..,-  -  Y  ese'  estiaates '  nave  tt*rr$-  obtained  for 'a 
rigid  Surface  and'  H'  :V#fore-  -m)tSs?5M'-yaluejs  that 

are.  1  ik*iy -'ro-'s>^eu.r,.  tir'-rh-v'  pressure  'itwirraara  on  n  deformable 
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surface,  the  localized  contact  zone  will  be  depressed;  hence, 
the  relative  velocity  of  approach  of  *  he  liquid  is  less  than 
what  it,  would  be  for  a  rigid  surface. 


Blowers  investigated  the  propagation  of 
stress-  waves  in  an  elastic  half-space  subjected  to  a  uniform 
pressure  P  distributed  over  an  expanding  circular  region  on 
'its  surface  whose  time-dependent  radius  is  deduced  from  the 
idealized  model  of  a  compressible  drop  striking  a  rigid  sur¬ 
face  depicted  in  Pig.  3^.  If  the  impact  velocity  is  Y_, 
then  the  radius  of  f  he  gon^ac*-  area,  a{*  ),  increases  with  time 
According  to  the  relation 


aft)  -v/?rv  ti-fV  t? 
o  o 


where  R  is  the  radius  of  the  spherical  drop.  The  radial  velo¬ 
city  of  the  expanding  contact  area  is 


( R- Vi  )V 

ypRvct-(v  t ) 


■w 


For  small  values  of  t  compared  to  the  duration  of  the  impact 


*  KV- 


The  pressure  distribution  on.  the  surface  of  ifte  solid  '{z-*?y) 
adopted  by  Flowers  ic  his  analysis  5s  of  to-  form 


P(r,t)  *  p  Tor  r<kt 
«  0  for  r>kt  ' 


Th'  magnitude  of  P  has  to  be  prescribed  and  is  taken  *o  be  a 
known  quantity. 


'  The  velocity  of  the  expanding  c< rcuiar  boundary 
is  initially  supersonic  but  decreases  as  t“\  The  disturbance 
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Figure  38.  Idealized  Model  fi  Evaluation  of  Time-dependf  nt  Radius  of  Expanding 
Contact  Area  for  a  Spherical  Drop  Striking  a  Solid  Surface 


experienced  by  ire-  solid  medium  will  then  be  due  to  the  pro¬ 
pagation  of  a  dilatation  wave  traveling  at  a  velocity  c^,  the 
propagation  of  a  distort j on al,  or  shear,  wave  traveling  at  a 
velocity  Cg,  and  :i  Rayleigh  wave  propagating  along  the  surface 
at  a  velocity  e.,  which  vr  slightly  less  than  c?. 

c  =  ^  E(!-r) 

1  /  p(  i+v)(l-?i»)  (39) 

„  _  rz e 

'2  /  2p7T+7)' 

where  E  is  Young's  modulus,  v  is  Poisson's  ratio,  and  pis 
the  density  uf  the  medium. 

As  long  as  the  boundary  of  the  loaded  region 
is  moving  at  a  velocity  A(t)^c1  the  radius  of  the  disturbance 
will  be  the  same  as  that  of  the  loaded  area.  When  S/tJcc^, 
the  disturbance  will  move  ahead  of  the  load  and  will  continue 
to  travel  at  the  sonic  velocity  of  the  medium.  The  time  at 
which  &('t)=c1  can  be  determined  from  Eq.  (37). 

T1  = 

X 

For  times  greater  than  r1?  Blowers  displayed  the  solutions  for 

the  development  of  the  dilatatl onal,  d istortional,  and  Rayleigh 

waves  as  they  propagate  ahead  of  the  loaded  area  along  the 

plane  surface  of  the  elastic  half-space.  Limited  numerical 

o 

results  are  provided  for  the  case  when  (cg/c^)1-  -  2/7  or 

0. 3.  This  value  of  Poisson's  ratio  is  approximately  equal 
to  thab  of  steel  for  which  -  19,500  fps.  Blowers  shows 
that  the  transient  radial  and  azimuthal  surface  stresses  at 
time  t  =  8^  are  primarily  compressive  with  extremely  high 
tensile  stresses  found  only  in  a  narrow  band  immediately  behind 
the  Rayleigh  wave  front.  Blowers  comments  that  although 
extreme  radial  tensions  arise,  the  extent  of  the  radial  dimen¬ 
sions  in  which  they  occur  is  always  quite  small.  Thus  each 
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point  experiences  the  critical  tension  for  an  extremely  short 
duration.  Blowers  conjectures,  that  this  result  of  the  analy¬ 
sis  tends  to  raise  doubt  in  the  effectiveness  of  the  Rayleigh 
wave  as  the  cause  of  liquid  impact  erosion  damage. 

Blowers'  analysis  assumes  that  the  magnitude 
of  the  applied  pressure  is  constant  over  an  expanding  circular 
area.  In  actuality  the  magnitude  of  the  maximum  value  in  the 
pressure  disiributton  and  the  form  of  the  distribution  is 
changing  with  time.  As  indicated  earlier  a  fluid  mechanics 
analysis  of  a  spherical  drop  striking  a  rigid  surface  has  been 
provided  by  Huang using  a  numerical  scheme.  The  results 
from  these  computations  agree  favorably  with  the  existing 
experimental  data.  Ine  pressure  distributions  at  the  liquid- 
solid  interface,  the  isobaric  distribution  within  the  drop, 
and  the  radial  velocity  of  the  contact  zone  have  been  computed 
as  a  function  of  a  non-dimensional  time  unit  and  the  spatial 
coordinates.  Computed  results  are  available  in  Huang's 
thesis^29)  for  a  2mm  wa+er  drop  impacting  at  9^0  fps.  Under 
these  conditions  one  dimensionless  time  unit,  =  1,  is  equal 
to  lgs.  For  these  impact  conditions  and  taking  the  .solid 
medium  to  be  steel  in  conformance  with  Blowers  computations, 
Eq.(4o)  yields 


=  4.24  x  10~9  sec. 
for  k  -  2.54  ft/(sec)^ 

The  time-dependent  radius  of  the  loaded  area 
computed  by  Huang  for  the  non- slip  boundary  condition  agrees 
exactly  with  the  radius  obtained  from  the  simple  model  used 
-lowers  and  given  by  the  relation  in  Eq.(37).  The  radii 
corresponding  to  the  full-slip  boundary  condition  also  employed 
by  Huang  differ  by  less  than  a  few  percent  from  the  values 
based  on  the  simple  model.  This  agreement,  is  found  to  hold  for 
values  of  the  contact  area  radius  approaching  the  radius  of 
the  drop. 
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The  computer  models  of  a  compressible  liquid 
drop  striking  a  rigid  surface  provide  a  reasonably  accurate 
approximation  to  the  pressure  distribution,  applied  to  a  rigid 
surface,  while  Blowers'  results  provide  a  quantitative  evalu¬ 
ation  of  the  transient  stress  distribution  within  an  elastic 
solid  subjected  to  a  known  uniform  pressure  distribution. 
Although  the  exact  conditions  prevailing  at  the  interface 
between  the  drop  and  deformable  solid  cannot  be  analyzed 
directly,  the  two  approaches  above  can  be  combined  in  a  way 
which  provides  at  least  a  reasonable  approximation  to  the 
coupled  problem.  For  extremely  low  impact  velocities  the 
pressure  distributions  obtained  from  the  droplet  models  can 
be  used  directly  to  ascertain  the  response  of  the  solid 
medium.  However,  at  impact,  velocities  which  are  sufficient 
to  compress  the  solid  by  a  significant  amount,  an  elementary 
shock  wave  analysis  of  the  interfacial  region  can  be  employed 
to  estimate  the  magnitude  of  the  pressure  transmitted  to . the 
solids 

In  reviewing  Huang's  results  it  is  seen  that 
the  pressure  distribution  is  hemispherical  up  to  the  point  at 
which  the  maximum  pressure  is  reached  and  then  it  assumes  a 
more  uniform  distribution  as  the  magnitude  of  the  pressure 
subsides  with  the  onset  of  substantial  lateral  outflow.  In 
order  to  employ  Blowers'  analysis,  the  applied  pressure  dis¬ 
tribution  during  both  periods  will  be  assumed  to  be  uniform 
and  have  a  magnitude  equal  to  the  average  value  of  the  actual 
distribution..  The  fluid  particle  velocities  corresponding  to 
this  pressure  can  be  determined  from  the  Hugoniot  relation 
for  water  at  each  instant  of  time.  The  procedure  is  illus¬ 
trated  in  Fig.  39.  A  family  of  curves,  reflection  Hugoniots 
for  water,  are  generated  and  plotted  as  shown  in  Fig.  39(b) 
with  time  as  a  parameter.  This  approach  allows  the  pressure 
pulse  to  have  a  finite  rise  time  in  conformity  with  the  two- 
dimensional  analysis  for  a  spherical  drop.  The  conditions 
prevailing  at  the  interface  between  the  drop  and  the  solid 


(b)  Time-dependent  Reflection  Hugoniots  for 
Water  with  Average  Time-dependent  Pressures 
Computed  from  Huang's  Analysis 


(c)  Time-dependent  Pressure  Applied  to  Solid 
Surface  Corresponds  to  Intersection  Points 
between  Hugoniot  for  the  Solid  and  the 
Reflection  Hugoniots  for  Water 


Figure  Tl>.  (-.valuation  of  Pressure  Transmitted  to 
a  Solid  Medium  by  l.i<|tiid  Droplet  Impacts 
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medium  are  approximated  by  plotting  the  pressure-particle 
velocity  relations  for  the  impacted  solid  on  Fig.  39(b)  as 
shown  in  Fig.  39(c).  Assuming  the  pressure  and  particle 
velocities  are  continuous  across  the  interface,  the  points  of 
interaction  of  the  reflection  Hugoniots  for  water  and  the 
Hugoniot  for  the  solid  provide  the  instantaneous  magnitude 
of  the  pressure  applied  to  the  solid  surface  accounting  for 
the  compressibility  of  the  solid  medium. 

A  fair  approximation  to  the  spatial-  and 
time-dependence  of  the  localized  pressure  imparted  to  a  deform¬ 
able  solid  by  a  spherical  liquid  drop  impacting  at  low  to 
moderate  velocities  has  now  been  achieved.  The  magnitude  of 
the  pressure  distribution  can  be  used  with  Blowers'  analysis 
to  evaluate  the  transient  stresses  which  arise  in  response 
to  the  disturbances  propagating  into  the  interior  of  an 
elastic  half-space.  As  the  magnitude  of  the  pressure  increases 
it  is  possible  that  finite  amplitude  waves  will  propagate  into 
the  solid  medium!33"35^ 

Recently,  Peterson ^3^  advocated  the  reduc¬ 
tion  of  the  transient  liquid  droplet  problem  to  a  quasistatic 
analysis  based  on  the  Hertzian  theory  of  impact  In  order  to 
evaluate  the  response  of  the  target  material.  On  the  basis 
of  Huang's  computer  analyses,  he  assumed  a  hemispherical  dis¬ 
tribution  of  pressure  as  representative  of  the  loading 
conditions  during  the  period  of  pressure  build-up.  Peterson's 
approach  would  be  valid  if  it  could  be  established  that  the 
quasistatic  stress  states  are  attained  quite  rapidly  after 
the  passage  of  the  Rayleigh  surface  waves. 

Once  again  we  refer  to  the  published  results 
of  Blowers  applicable  to  a  uniform  load  applied  to  an  expand¬ 
ing  contact  area  representative  of  a  spherical  droplet  colli¬ 
sion.  For  a  very  short  time  interval  after  the  Rayleigh  wave 
propagates  ahead  of  the  moving  periphery  of  the  loaded  area, 
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both  the  radial  and  circumferential  stress  components  are 
tensile  and  of  unequal  magnitude  in  the  region  between  the 
Rayleigh  wave  front  and  the  boundary  of  the  contact  area  on 
the  surface  of  the  naif  space.  The  static  solution  for  a 
uniform  load  acting  over  a  circular  region  on  the  surface  of 
an  elastic  half-space  shows  that  the  magnitude  of  the  radial 
stress  and  circumferential  stress  components  are  equal  but 
the  circumferential  stress  is  compressive  while  the  radial 
stress  is  tensile  in  the  region  exterior  to  the  loaded 
area. (37)  jn  principle  it  is  possible  to  evaluate  the  form 
of  the  stress  components  at  larger  time  intervals  using 
Blowers'  results,  however  access  to  his  computer  program  has 
not  been  established  at  the  present  time.  The  validity  of 
Peterson's  hypothesis  t.ha<  the  static  stress  distributions 
can  be  used  after  the  passage  of  the  Rayleigh  wave  to  esti¬ 
mate  the  response  of  *he  target  materials  to  droplet  impacts 
can  therefore  be  tested  directly.  Both  Blowers and 
Eason point  out  that  the  static  values  of  the  displace¬ 
ments  are  approached  quite  rapidly  along  the  axis  of  symmetry 
of  the  loaded  area. 

A  significant  amount  of  additional  work  is 
required  to  obtain  an  operational  computer  code  which  will 
provide  meaningful  quantitative  results.  However,  such  a  pro¬ 
gram  would  properly  delineate  the  conditions  under  which  a 
static  analysis  is  acceptable  and  those  for  which  a  transient 
stress  wave  problem  must  be  considered.  The  three-dimensional 
elastic  analysis  also  provides  more  accurate  evaluations  of 
the  magnitude  of  the  stress  components  within  the  solid  medium 
than  the  elementary  one-dimensional  approaches  that  have  been 
prevalent  in  erosion  studies.  Having  more  accurate  approxi¬ 
mations  to  the  maximum  stress  state'  the  solid  must  endure  and 
their  duration  will  lead  effectively  to  the  development 

of  realistic  connections  with  a  matt  rial's  mechanical  and 
physical  properties,  fracture  behavior  and  microstructurc. 
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3.5.2  Experimental  Erosion  Data 

vJhen  the  same  specimen  holder  used  for  the 
solid  particle  impacts  was  used  in  the  rain  erosion  tests 
without  additional  modification,  the  glass  specimens  exhibited 
a  strong  tendency  for  cracking  in  the  early  stage  of  each  test 
A  number  of  changes  were  made,  however  only  a  limited  amount 
of  experimental  data  pertaining  to  liquid  drop  impacts  were 
obtained  within  the  present  program.  In  subsequent  testing 
smaller  specimens  wilt  be  employed  set  in  a  thick  polyurethane 
support  Which  is  inserted  into  the  existing  specimen  holder 
at  th-  end  of  the  rotating  arm.  No  weight-loss  data  were 
obtained,  since  the  specimen  fractured  before  more  than  a  few 
weight  measurements  could  be  made. 

Tne  cumulative  damage  on  the  surface  of 
specimen  PY-1  was  negligible  after  an  exposure  time  of  40 
min.  An  occasional  microscopic  crack  could  be  found  on  the 
surface,  but  it  appears  that  the  impact  velocity  of  300  fps 
is  not  sufficient  to  produce  extensive  surface  cracking. 

The  fine  cracks  which  did  appear  did  not  show  any  particular 
characteristic  which  could  be  related  to  the  pressure  dis¬ 
tribution  generally  assumed  to  be  associated  with  water  drop¬ 
let  collisions.  A  small  number  of  flaws  must  be  present  on 
the  surface  which  are  of  sufficient  size  to  extend  to  low- 
power  microscopically  observable  cracks  at  t tv  pressures 
developed  by  a  drop  Impacting  at  300  fps. 

Specimen  FY-2  showed  a  sparse  distribution 
of  microscopically  visible  cracks  penetrating  the  interior 
of  the  specimen.  There  was  some  tendency  toward  a  concentric 
ring  configuration  with  at  least  a  portion  of  the  cracks 
centered  around  an  identifiable  origin.  The  slope  of  the 
fracture-  surfaces  extending  into  the  specimen  provided  some 
idea  of  the  direction  of  the  transient  stress  pulses  generated 
by  droplet  collisions.  <"io3ed  curves  where  the  internal 
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crack  Intersects  the  surface  were  found  on  the  surfaces  of 
specimens  PY-2,  PY-3  and  PY-5*  The  curves  were  quite  irregular 
and  did  not  have  the  circular  form  observed  by  Bowden  and 
Field for  jets  striking  glass  surfaces  at  considerably 
higher  impact  velocities.  The  lateral  dimensions  of  the  closed 
regions  ranged  from  a  few  mils  up  to  13  mils.  Some  areas 
did  show  a  slight  tendency  toward  concentric  ring  formation 
within  pie-shaped  regions  with  an  included  angle  in  the  range 
of  45  to  60  degrees.  Several  cracked  regions  on  specimen  PY-3 
originated  at  semi-circular  cracks  with  radii  around  2  to  6 
mils.  The  cracks  penetrated  60  to  80  microns  into  the  bulk 
material.  .It  is  estimated  that  the  fracture  surfaces  are 
inclined  at  angles  from  29  to  38  degrees  to  the  surface  of 
the  specimen. 

The  majority  of  the  subsurface  fractures  had 
a  characteristic  form.  The  fracture  surface  extending  from 
the  surface  of  the  specimen  into  the  bulk  material  was  very 
smooth  for  a  sh eft  distance  below  the  surface  and  then  exhi¬ 
bited  a  hackled  appearance  for  the  remainder  of  its  extent. 

Fine-scale  chipping  initiates  through  the 
intersection  of  the  random  crack  formations  on  the  surface. 
Specimen  PY-5  gave  some  indication  of  the  longer-term  erosion 
behavior.  Four  large  pits  which  were  0.1  to  0.2  in,  wide 
were  present  on  the  surface  before  the  specimen  fractured. 

At  500  fps  the  material  removal  mechanisms  are  quite  pre¬ 
ferential  and  the  erosion  process  conforms  to  the  fundamental 
hypothesis  of  the  erosion  model. 
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4.0  GENERAL  EROSION  MECHANISMS 

Table  X  lists  the  erosion  tests  conducted  on  ultra-high 
molecular-weight  polyethylene,  a  nickel-cnated  epoxy  composite 
and  a  graphite-fiber  reinforced  epoxy.  The  specimens  were 
exposed  to  a  standard  rainfield  of’  1  in/hr  of  1.8  mm  drops 
in  the  AFML- fell  erosion  facility.  The  90-micron  sand  des¬ 
cribed  in  Section  3*1  was  used  for  the  solid  particle  erosive 
mdeium.  In  all  nf  the  >ests  the  flow  rate  was  260  g/'min. 

4.1  Erosion  of  Ultra-High  Molecular-Weight  Polyethylene 

A  series  of  flat-faced  and  standard  aerodynamicaliy- 
shaped  polyethylene  specimens  were  subjected  to  the  standard 
rainfield  (1  in /hr.  of  1.8  mm  droplets)  and  the  90-micron  sand 
characterised  in  Fig.  2  and  3*  The  ultra-high  molecular  weight 
(UHMW)  polyethylene  is  Host algen  GUR  supplied  by  the  American 
Hoechst  Corporation.  The  flat -faced  polyethylene  specimens 
were  rectangular  parallelepipeds  3*25  in.  x  0.5  in. 

The  magnitude  of  the  eroding  area  is  different  on 
the  surface  of  each  specimen  configuration.  The  eroding  area 
for  the  flat- faced  specimens  is  2,25  sq.  in,  (3.0  in.  x 
0.75  in. }  and  for  the  ae*odynamic  shape  it  is  1.25  sq.  in. 
However  there  is  an  angular  effect  due  to  the  curvature  of 
the  aerodynamic  shape  which  makes  a  direct  comparison  between 
each  specimen  configuration  somewhat  difficult.  The  flat-faced 
specimens  were  used  in  an  attempt  to  eliminate  the  influence 
of  this  curvature  on  the  erosion  rates  and  to  facilitate  the 
microscopic  examination  of  the  e  rcsion  damage. 

The  surfaces  of  the  polyethylene  specimens  were 
polished  before  testing.  The  condition  of  ‘he  initial  sur¬ 
faces  before  testing  IS  shown  in  Fig.  40.  A  point  was  reached 
in  the  polishing  operation  beyond  which  additional  polishing 
of  the  surface  did  not  produce  any  refinement  *he  surface 
condition  on  a  microscopic  level. 
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■  TABLE  X 

EROSION  TEST  SCHEDULE  -  OENERAL  EROSION  MECHANISM  STUDIES 
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4,1.1  .Rain  F "c-p ion 

Rain  erosion  tes's  w  three  specimens  of 
polyethylene  at  500,  730,  and  1150  fps  indicated  that  the 
time  scale  for  which  . a  particular  erosion  mechanism  is  opera¬ 
tive  ih  thebe  specimens  is  only  contracted  or  expanded  an  the 
impact  velocity  is  increased  or  decreased  but  no  change  is 
evident  in  the  ■  sequence  f  operat ive  erosion  mechanisms .  It 
was  therefore  decided  to  conduct,  the  remaining  controlled 
erosion  tests  at  a  single  velocity  instead  of  a  range  of  velo¬ 
cities  in  order  to  concentrate  on  describing  the  erosive 
•bfhaVlor  in  more  detail  ♦ban  would  be  possible  if  each,  -s.pecir 
men, was  run  a*  a  different  velocity;  "his  approach  is  adopted 
in  all  of  the  investigations,  of  •  erosion.  mechanisms  for.  the  ■  ' 
material:  systems  considered  in  this  research,  As  .long  as  the 
form  of  the  erosion,  mechanisms  (operative  in  these  material 
■systems ^  are  not  dependent  ort  t  he  particle  impact •■•'yeloc.it  jU*s.. 
the  vs?  loci  ty- depends  no  o  if  the  eroslcn  rates  can  be  easily.,;  ;.•• 
.prescribed  after  the  sequroeg  vf  *»ros ion  mechanisms  in  a 
piS*t icular  material  system., have  bwn; Identified  and  described. 
ara.lyUca.Uy', 

.  .  for  the  am'dynamleally- shaped  specimens 
exposed (to.  the  #'t an&ard ■Tiinne.id  at  impact  velocities  of  '-'*00., 
:i ..and  1*20  fps  it  wa^  fovtnd  . that  at  500  fps  nc  weight  loes 
was  recorded  an  exposure  time  of  30  .minutes  (specimen 

?$)  and  at  1120  fps  the  weight  loss  was  too  large  after  30 
seconds  for  detailed  investigations  of  the  evolving  erosion 
mechanisms.  Therefore  73$  fps  war  selected  as  ?  he  moat  ?n t  *  - 
aMe  velocity  for  obtaining  weight -loss  data  and  for  inv««*  i- 
gat irg  the  erosion  mechanisms  responsible  for  material  removal 
At  500  fps  no  weight  loss  was  »ecorded  for  thv  flat-faced 
specimen,  '?*»  after  a  60-minute  exposure  time,  however  &  low 
density  of  small  pits,  less  than  5  mils  it.  diameter, 
beg  inning  to  develop  on  the  surface  at  “Ms  Mae. 
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Weight-loss  curves  were  obtained  for  bott: 
specimen  configurations  at  73°  fps  as  shown  in  Fig.  4l.  The 
data  are  net  corrected  for  differences  in  the  areas  exposed 
to  the  erosive  environments. 

At  730  fps  a  uniform  distribution  of  fibers 
were  raised  from  the  surface  of  specimens  e8  to  E10  after  a 
very  short  exposure  to  the  rain  environment.  After  an  expo¬ 
sure  of  three  minutes  small  pits  began  to  appear  on  the  sur¬ 
face  of  specimens  E8  and  E9.  The  average  size  of  the  pits  i  s 
five  mils.  The  mechanism  by  which  the  pits  are  formed  appears 
to  be  the  raising  of  portions  of  the  material  surrounding  a 
small  crack  in  the  material.  Several  pits  are  devoid  of  this 
raised  portion,  but  it  has  apparently  been  torn  away  during 
subsequent  impacts.  The  inferior 'of  the  small  pits  show  a 
fibrous' composition  along  the  walls  hut  the  base  of  ►he  pit 
is . generally  smooth.  The  pit  density  on  the  surface. of  -he 
specimen  is  not  very  great.  A  strong  tendency  for  the  .pita 
which  were  formed  to  grow  was  clearly  recognized  in  the  sub¬ 
sequent  micrographs.  Sequences  of  growth  in  these  specimens 
can  he  s-.en  in  Fig.  'IP  hh.  permanent  .deformations  on  the 
surface  cf  the  speeHms  were  also  observed  after  at;,  exposure 
time  of  ten  minutes.  Limited  weight-loss  data  have  been 
obtained  for  these  specimens.  The  surface  of  the  specimens 
.at  the  md  of  the  tests  were  very  **ougb.,  nowev».- r  It  did  not 
appear  that  any  new  mechanisms  were  being  activated  at  this 
time.  The  preferential  growth  of  pitted  regions  was  still 
prevalent.  Prior  to  this  *  ime  it  was  observed  that  the  line 
depressions  in  the  surface  did  provide  si‘es  for. ‘he  addi- 
t tonal  removal  of  material  by  tearing  pieces  from  ♦he 
surface.  This  erosion  process  'occurs  after  the  formation  of 
discrete  pits  and  ran  be  In  Pig.  99(e). 

9.1.2  hand  Erosion 

Use  sand  erosion  of  t ii<-  aerodynamic- shaped 
specimens  of  polyethylene  was  invest  1  -a*  «•.<*  using  the  uo-mJeron 
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(b)  Exposure  time  15  min.  (110X) 

Figure  i<3  Long-time  pit  growth  on  the  surface  of  polyethylene 
specimen  E9  (730  fps  in  standard  rainfield). 


P 

*  *-  J  * 


(d)  Exposure  time  40  min.  (60X ) 


Figure  43 


Long-time  pit  growth  on  the  .surface  of  polyethylene 
specimen  E9  (730  fpsin  standard  rainfield). 
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(a)  Exposure  time  7  min.  (110X) 


Figure  M  Pit  growth  and  coalescence  on  the  surface  of 
polyethylene  specimen  E9  (730  fps  in  standard 
rainfield). 


Figure  Uj  pit  growth  and  coalescence  on  the  surface  cf 
polyethylene  specimen  KQ  (730  fpr,  in  standard 
rainf ield ) . 


sand  at  a  flow  rate  of  260  g/mln.  and  an  impact  velocity  of 
880  fps.  The  weight-loss  data  are  plotted  in  Fig,  45.  It 
is  seen  that  a  steady  rate  of  erosion  is  quickly  established. 
This  data  can  be  compared  with  that  previously  reported. (39) 


The  sand  erosion  mechanism  is  quite  distinct 
from  that  associated  with  rain  droplets.  The  permanent 
indentations  of  the  sand  particles  are  readily  identified  along 
the  leading  edge  of  the  specimen  and  the  cutting  action  of  the 
particles  can  be  seen  along  the  curved  sides  of  the  specimen. 
Most  of  the  material  removed  from  the  specimen  occurs  along 
the  sides  of  the  specimen  between  the  angles  of  20°  and  50° 
as  measured  from  a  plane  passing  through  the  longitudinal 
axis  and  leading  edge  of  the  specimen. 


4.2  Erosion  of  Nickel-Coated,  Glass-Fiber  Reinforced 
Epoxy 

A  series  of  ? lx  specimens  of  a  nickel-coated,  glass- 
fiber  reinforced  epoxy  with  coating  thicknesses  of  3  and  6  mils 
were  exposed  to  rain  and  sand  erosion  environments.  The  pre¬ 
paration  of  the  specimens  is  described  by  Weaver  The 

specimens  tested,  as  recorded  in  Table  X,  are  the  soft-nickel 
coatings, 


The  initial  cond  it i on  of  the  surface  of  these 
specimens  is  indicated  in  Fig.  46.  The  surfaces  are  moder¬ 
ately  rough  with  *he  6  mil  coatings  being  coarser  than  the 


3  mil  coatings.  Two  types  of  initial  surface  features  are 
observed:  small  clyindrical  protrusions  and  extended  pits. 


4.2.1  Rain  Erosion 

Two  specimens,  E-l-3  with  a  3-mil  coating 
and  D-4-6  with  a  6-mil  coating,  were  eroded  at  730  fps  in 
the  standard  ralnfield,  it  was  evident  from  optical  microsopie 


F r.tWRE  Weight  lose  for  direct  impacts  of  90-rmicror,  sand  on 

polyethylene  at  880  fps. 


(a) 


Specimen  ..P-3- 


3 


(b)  Specimen  D- ?~3 


Figure  ;*£  Initial  surface  condition  of  nickel-coated 

glass  fiber  reinforced  epoxy  specimens. (110X ) 


examination  that  t^e  surfaces  were  experiencing  very  r 3 ne- scale 
plastic  deformation  without  pitting.  The  initial  phase  of 
the  erosion  process  produced  a  distribution  of  cracks  on  the 
surface  which  only  extend  a  portion  of  the  way  into  the  thick¬ 
ness  of  the  coating.  These  fractures  generally  occur  most 
readily  at  the  base  of  the  protrusions  which  are  initially 
present  on  the  s  urface.  The  cracking  becomes  more  extensive 
as  the  exposure  time  increases.  Over  the  duration  of  these 
tests  no  weight  loss  was  recorded.  On  a  microscopic  level 
there  was  no  evidence  of  any  erosion  mechanisms  which  would 
contribute  to  material  removal. 

Failure  occurred  in  specimen  P-1-3  at  an 
exposure  time  of  15  minutes  by  localized  buckling  of  the 
coating  at  the  outboard  end  of  the  specimen.  The  overall  con¬ 
figuration  of  the  failure  is  shown  in  Fig.  k?»  The  buckled 
ridges  extend  300  tc  **00  microns  above  the  surface  of  the  coat¬ 
ing.  Fig.  k8  and  k9  show  the  cracked  region  surrounding  the 
ridges  as  well  as  the  development  of  pits  which  extend  well 
Into  the  coating  and  in  some  cases  to  the  surface  of  the  sub¬ 
strate. 

In  order  to  se*  what  would  happen  to  the 
coating  now  that  it  had  buckled  locally,  the  apecume;. 
placed  on  the  rotating  arm  and  monitored  for  evide^e;.^  ^ross 
damage.  The  test  was  terminated  at  19  aihutes  wh*nO£.j\crge 
portion  of  the  coating  surrounding  the  buckled  re  sic:/  we* 
completely  removed  from  the  substrate  as  shown  in  Fig.  50. 

At  this  time  seme  pitting  and  extensive  cracking  was  observed 
near  the  center  or  the  specimen  as  illustrated  in  Hg.  51, 

Specimen  B-k- 6  with  a  6-mil  costing  was 
found  to  be  considerably  more  resistant  to  the  erosive  environ** 
went,  negligible  damage  war  observed  at  a  microscopic  level 
until  an  exposure  time  of  kO  minutes  when  cracking  began  to 
develop  on  the  surface.  The  severity  of  the  cracking  increased 
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Figure  49  Condition  of  the  surface  near  the  center  of 
specimen  D-l-3  after  15-minute  exposure  to 
the  standard  rainf.teld  at  730  fps,  (110X) 


-1 


with  exposure  time  and  a  low  density  of  deep  pits  were  begin- 
ning  to  evolve.  Microscopic  examination  of  the  specimen 
showed  that  regions  with  lateral  dimensions  of  10  to  20  mils 
were  being  depressed  or  raised  above  the  original  surface  at 
the  inboard  end  of  the  specimen  at  an  exposure  time  of  100 
minutes.  Some  idea  of  the  surface  conditions  at  this  time 
can  be  obtained  from  Pig.  52  for  the  inboard  end  and  from 
Fig.  53  near  the  center  of  the  specimen.  The  specimen  was 
carefully  monitored  as  the  time  of  exposure  was  further  in¬ 
creased,  until  a  large  hole  developed  in  the  coating  at  an 
exposure  time  of  108  minutes  as  shown  in  Fig.  54.  The  test 
was  then  terminated/ 

Although  the'  findings  made  on  the  basis  of 
these  two  specimens  will  have  to  be  confirmed  on  the  basis 
of  additional  tests,  it  may  be  conjectured  from  these  tests 
that  the  mode  of  failure  depends  on  the  thickness  of  the  coat¬ 
ing.  The  droplet  impacts  do  not  appear  to  produce  localized  . 
pitting  on  the  surface  as  is  the  case  of  bulk  nickel  specimens;, 
they  tend  rather  to  deform  the  coating  plastically  as  seen  by 
the  shallow  depressions  formed  on  the  surface  and  slowly  weaken 
the  bond  at  the  interface  between  the  coating  and  the  sub¬ 
strate.  The  3-roil  coating  o_nce  freed  from  the  substrate  does 
not. possess  enough  rigidity  to  be s 1st  the  centrifugal  forces 
applied  to  the  specimen  in  the  rotating  arm  appratus  and  so' 
buckling  of  the  surface  coating  results.  The  6-mil  coating 
fails  in  a  different  way,  since  it  is  capable  of  resisting  to 
a  greater  degree  the  applied  centrifugal  loads.:  The  rigidity 
of  the  coating  and  the  applied  centrifugal  forces  are  not 
usually  a  consideration,  however  the  ‘"3-mil  coating  must  not 
be  thick  enough  to  resist  the  loads  imposed  by  the  erosion 
apparatus  even. at  low  velocities. 

4.2.2  Sand  Erosion 

The  failure  mechanism  for  the  nickel-coated 
specimens  eroded  by  the  90-micron  sand  was  quite  different 
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Figure  52  Condition  of  the  coating  surface  at  the  inboard 
end  of  specimen  D-4-6  after  100-minute  exposure 
to  the  standard  rainfield  at  730  fps.  (110X) 


Figure  53 


Condition  of  the  coating  surface  at  the  center  of 
specimen  D-4-6  after  100-minute  exposure  to  the 
standard  rainfield  at  730  fps.  (15X) 


Figure  5^  Failure  of  the  coating  at  the  inboard  end  of 

specimen  D-4-6  after  108-minute  exposure  to  the 
standard  rainfield  at  730  fps.  ( 15X ) 
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from  that  found  for  the  standard  rainfield.  Two  specimens 
with  a  3-mil  coaiin-  (P-7-3  and  D-3-3^  and  two  specimens  with 
a  6-mil  coating  ( D- - 1  and  D-6-6)  were  exposed  to  sand  erosion 
at  730  fps  as  indicated  in  Table  X.  In  this  cas--  there  was 
a  measurable  w-'ight  loss  from  the  coating  surfaces  as  summar¬ 
ized  in  Fig.  55  and  56.  It  is  seen  from  these  figures  that 
after  a  relatively  short  incubation  period  a  uniform  rate  of 
material  removal  is  established.  The  sand  flow  rate  was  260g/ 
min . 


For  both  coating  thicknesses  failure  occurred 
by  buckling  of  the  coating  which  was  initiated  at  the  out¬ 
board  end  of  the  specimen.  No  cracking  of  the  surface  was 
observed  as  in  the  case  of  rain  erosion;  however,  initiation 
of  buckling  of  the  coating  took  place  at  an  exposure  time  of 
i<5  seconds  for  the  3-mil  coatings  and  approximately  85  seconds 
for  the  6-mil  coatings.  After  the  initial  buckling  the  condi¬ 
tion  of  the  surfaces  progressively  worsened,  until  they  reached 
the  conditions  shown  in  Fig.  57  and  58  at  the  end  of  each  test. 

It  is  reasonable  to  assume  that  the  reduction 
in  the  coating  thickness  is  critical  for  failure  in  the  case 
of  sand  erosion  and  that  buckling  of  the  coating  can  only  occur 
if  it  is  separated  locally  from  the  substrate.  Specimens 
D-2-3  and  D-6-6  were  eroded  by  the  90-micron  sand  for  a  total 
exposure  time  slightly  less  than  that  at  which  initial  buckling 
of  the  coating  first  occurred  in  specimens  D-3-3  and  D-5-6, 
respectively.  These  specimens  were  cut  with  a  diamond  wheel 
at  a  very  low  rate  of  penetration,  so  the  bonding  between  the 
coating  and  substrate  would  not  be  disturbed.  The  sections 
were  mounted  and  polished.  Microscopic  examination  of  the 
interface  in  all  cases  revealed  that  the  integrity  of  the 
bond  was  unaffected  by  the  eroding  forces.  No  evidence  for 
the  initiation  of  any  kind  of  failure  could  be  seen. 
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:-3-3  eroded  by  90-mic 


FIGURE  56  Weight  loss  for  specimen  D-5-6  eroded  by  90-micron  quartz 
sand  at  730  fps. 


(a)  1.75X  I 

f 
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'J8  Buckled  surface  of  specimen  D-5-6. 
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The  interface  for  cross  sections  of  segments 
from  specimens  D-3-3  and  D-5-6  away  from  the  buckled  region 
were  also  examined  microscopically  and  found  to  be  secure. 

The  mechanism,  or  critical  conditions,  for  failure  of  the 
coated  composites  has  not  yet  been  discovered. 


4.3  Erosion  of  Graphite-Fiber  Reinforced  Epoxy 

Aerodynamic- shaped  specimens  were  exposed  to  the 
standard  rainfield  and  the  90-micron  sand  at  a  flow  rate  of 
260  g/min.  HT-S  graphite  fibers  (manufactured  by  Hercules 
Incorporated)  were  placed  in  a  matrix  of  ERL  2256  epoxy 
(Union  Carbide  Corporation)  with  ZZLB-0820  hardener.  The 
material  properties  of  these  components  of  the  composite  are 
provided  in  Table  XI.  The  specimens  were  obtained  from  the 
Dayton  Research  Institute,  Composite  and  Polymeric  Materials 
Laboratory. 

The  method  of  preparation  was  as  follows: 

The  epoxy  resin  was  prepared  by  mixing  27  parts  of 
ZZLB-0820  hardener  per  100  parts  of  ERL  2258  epoxy.  This 
mixture  was  diluted  with  75  parts  of  methyl  ethyl  ketone  per 
100  parts  cf  epoxy  to  facilitate  impregnating  the  fiber  on 
a  filament  wii  ding  machine.  The  epoxy  impregnated  fiber  was 
cut  to  the  mold  size  and  48  plies  of  material  were  stacked  in 
the  mold.  Curing  of  the  specimens  was  accomplished  by  the 
following  schedule:  2  hr.  at  l80°F  followed  by  4  hr.  at  300°F 
in  a  laminating  press.  The  mold  was  then  cooled  to  ambient 
temperature  under  pressure. 

Microscopic  examination  of  the  original  surfaces  of 
specimens  indicated  they  were  smooth  except,  for  numerous 
depressions  where  an  area  of  the  graphite  fibers  were  missing 
due  to  fabrication  of  the  composite.  The  condition  of  the 
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TABLE  XI 


MATERIAL  PROPERTIES  OF  COMPONENTS  OF 
GRAPHITE- FIBER  REINFORCED  EPOXY  COMPOSITES 


Tensile  strength  (psi) 
Tensile  modulus  (psi) 
Tensile  elongation  (%) 

Heat  distortion  temp.  (°C) 
Density  (g/cc) 


ERL2256/ZZLB  0820 

( 27  pph ) 

1.4  x  10^ 

4.2  x  105 

6.5 

145  -  150 


HT-S 

Graphite 

Fibers 

3.71  x  105 

41.06  x  106 


1.8 
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original  surface  is  shown  in  Fig.  59*  The  appearance  of  the 
surface  indicates  that  a  moderate  porosity  may  prevail  within 
the  interior  of  these  specimens.  In  order  to  determine  the 
internal  constitution  of  the  graphite-fiber  reinforced  epoxy 
specimens,  the  four  specimens  were  sectioned  and  examined 
microscopically.  It  was  found  that  a  high  porosity  extended 
throughout  the  interior  of  the  specimens  as  can  be  seen  in 
the  tynical  cross  sections  in  Fig.  60.  The  porosity  was  esti¬ 
mated  to  be  on  the  order  of  five  percent.  The  fiber  content 
for  these  composites  was  determined,  to  be  31  percent  by  volume. 
The  fibers  are  distributed  nonuniformly  throughout  the  com¬ 
posite  as  clearly  seen  in  Fig.  60(c)  and  (d). 

The  particle  impact  velocity  was  held  constant 
at  730  fps  in  the  erosion  tests.  The  duration  of  each  test 
and  its  objective  are  listed  in  Table  X. 

4.3.1  Rain  Erosion 

Thf  erosion  characteristics  of  specimen  C-l 
exposed  to  the  standard  rainfield  are  quite  distinct  from  those 
for  solid  particle  erosion.  Considerably  larger  fragments 
of  material  are  removed  in  this  case  due  to  cracking  of  the 
matrix  material.  The  general  features  of  various  areas  on  the 
eroded  surface  of  specimen  C-l  are  shown  in  Fig.  6l  to  64. 
Details  of  the  fractured  fibers  can  be  seen  in  Fig.  6l(b) 
and  62(b).  Thr-  appearance  of  the  eroded  surfaces  indicates 
that  bundles  of  fibers  and  the  surrounding  matrix  are  removed 
as  short  longitudinal  segments  as  observed  in  Fig.  6l(a).  The 
missing  segment  above  the  deep  trough  in  the  center  of  the 
micrograph  which  is  about  10  mils  long  and  three  fiber  dia¬ 
meters  thick  with  the  adjacent  segments  which  have  also  been 
removed  roughly  half  this  length.  Cracks  in  the  matrix  mate¬ 
rial  are  clearly  visible  in  all  these  micrographs  and  a  broad, 
deep  crack  extending  under  a  larger  volume  of  the  composite 
material  is  shown  in  Fig.  64.  Evidence  for  the  removal  of  a 
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short  segment  -approximately  6  mils  in  length  is  shown  in  the 
center  of  Fig.  62(a)  with  a  detail  of  the  fracture  surface  pro¬ 
vided  in  Fig.  62(b). 


A  preliminary  hypothesis  for  the  operative 
erosion  mechanism  for  the  intermediate  to  later  stages  of 
the  erosion  process  has  been  formulated  on  the  basis  of  these 
microscopic  investigations.  Thb  direct  impact  of  the  water 
droplets  produce  local  fractures  which  penetrate  the  matrix 
material  but  do  not  crack  the  graphite  fibers.  Separation 
occurs  at  the  fiber/matrix  interface  which  would  indicate  the 
bonding  between  the  fiber  and  matrix  is  the  critical  component 
in  the  composite  structure.  Once  a  longitudinal  crack  has 
been  nucleated  parallel  to  the  direction  of  the  fibers  it  will 
grow  during  subsequent  impacts  in  addition  to  the  formation 
of  other  longitudinal  cracks  in  the  vicinity  of  the  first  one. 

A  longitudinal  fracture  network  will  continue  to  develop  until 
segments  of  the  composite  are  nearly  freed  from  the  bulk  of 
the  material.  The  pressure  from  a  properly  oriented  droplet 
applied  to  the.  cracks  surrounding  a  material  segment  will 
push  the  segment  out  from  the  surface  and  produce  intense 
localized  bending.  This  condition  should  be  roughly  equiva¬ 
lent  to  a  beam  rigidly  supported  at  its.  end  and  subjected  to 
a  transverse  distributed  load.  This  observation  provides  the 
basis  for  a  simple  model  which  may  be  developed  from  elemen¬ 
tary  concepts,  from  beam  theory.  The  micrographs  of  the 
fracture  planes  at  the  ends  of  these  segments  (Fig.  61(b)  and 
63(b)  show  that  the  fibers  behave  in  a  perfectly  brittle  manner. 
The  fracture  network  can.  extend  beyond  the  ends  of  the  initial, 
segment  removed  making  it  possible  for  removal  of  additional 
fragments  along  a  longitudinal  bundle,  of  fibers.  At  this 
point  the  local  conditions  can  be  considered  equivalent  to 
a  cantilever  beam  subjected  to  a  distributed  load.  This  pro¬ 
cess  can  continue  along  a  group  of  fibers  until  conditions 
become  unfavorable  to  further,  fragmentation  as  in  the  case  of 
the  large  fiber  bundle  consisting  of  seven  to  eight  layers 
in.  1  he  center  of  Fig.  61(a)'. 


The  micrographs  in  Fig.  65  show  voids  in 
specimen  C-l  which  are  beginning  to  be  uncovered.  The  frac¬ 
tures  at  the  ends  of  these  voids  are  similar  to  those  observed 
for  regions  in  which  no  voids  are  found  to  exist.  The  signi¬ 
ficance  -of  the  voids  in  relation  to  the  overall  erosion  pro¬ 
cess  is  not  known  at  this  time,  however  in  the  case  of  rain 
erosion  they  should  provide  favorable  conditions  for  subsur¬ 
face  cracking  of  the  epoxy  matrix. 

The  removal  of  material  from  the  surface  of 
the  specimens  by  both  fluid  and  solid  particle  impacts  takes 
place  by  the  removal  of  pieces  of  the  complete  composite  in 
a  manner  similar  to  a  brittle,  homogeneous  material.  The 
presence  of  the  fiber  reinforcement  does  not  improve  the  ero¬ 
sion  resistance  of  the  composite.  In  the  case  of  rain  erosion 
progressively  larger  pieces  of  material  are  removed  from  the 
specimen  as  the  exposure  time  increases.  The  bonding  surfaces 
of  the  fibers  lying  parallel  to  the  longitudinal  axis  of  the 
specimen  become  more  and  more  susceptible  to  penetration  and 
cracking  by  the  rain  drop  impacts  as  the  erosion  pits  in  the 
surface  deepen  with  the  length  of  exposure  to  the  erosive 
environment.  This  effect  is  clearly  seen  in  the  shape  of 
the  weight-loss  curve  in  Fig.  66.  The  condition  of  the  sur¬ 
face  of  specimen  C-l  at  the  end  of  the  weight-loss  evaluation 
is  shown  in  Fig.  67. . 

4,3.2  Sand  Erosion 

The  condition  of  specimen  C-4  when  exposed 
to  a  sand  environment  is  shown  in  the  micrographs  in  Fig.  68. 
Material  is  removed  on  a  very  fine  scale  by  localized  crack¬ 
ing  of  both  the  fibers  and  matrix.  The  exposed  fibers  are 
relatively  free  of  the  surrounding  epoxy  (Fig.  69),  so  t he 
epoxy  must  be  chipped  free  from  the  surface  of  the  fibers 
and  then  the  exposed  fibers  are  pulverised  during  subsequent 
impacts.  Aerodynamic- shaped  specimens  were  used,  so  the 
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erosion  damage  will  be  a  combination  of  direct  deformation 
and  wear  deformation.  The  micrographs  in  Pig.  68  and  69 
depict  the  surface  condition  near  the  leading  edge  of  the 
specimen,  while  the  micrographs  in  Fig.  70  show  the  surface 
condition  away  from  -lie  leading  edge  where  the  overlying  epoxy 
Is  still  being  removed.  The  micrographs  for  the  erosion  by 
..Qtumicron  sand  particles  on  specimen  C-k  in  Fig.  68(a)  shew 
areas  where  well-defined  particle  imprints  can  be  identified. 
The  fracture  boundary  crossing  a  number  uf  fibers  shows  a 
degree  of  uniformity  which  would  indicate  that  it  was  produced 
by  the  shearing  action -of  an  impacting,  irregular- shaped  sand 
particle.  The  nrimary  erosion  meci  anism'in  this  case  is  frag¬ 
mentation.  of  the  matrix  and  fibers  which  occurs  at  a  uniform 
rate  throughout  the  entire  exposure  time  except  for  the 
Initial  Mage  of  the  erosion  .'process  as  evidenced  by  the 
linearity  of  the  weight-loss  curve  In  Fig.  71.  The  smooth¬ 
ness  of  the  eroded  surface  after  100  mg  weight  loss  can  be 
seen  in  Fig.  72.  The  sand  erosion  of  specimens  c4.to  c-6 
did  not  exhibit  the  dependence  on  the  angle  of  attack  for  -the 
impacting  particles  which  was  found  to  prevail  for  the  poly¬ 
ethylene  specimens,  7h*  degree  of  material,  removal  war  con¬ 
stant  up  to  an  angle  of  6ptr  as  measured  from  'he  plane  passing 
through  the  longitudinal  axis  and  t.n«  leading  edge  of 
aeredynami c- shaped  spec!  *cn, 

-The  micrograph's  in  Fig.  7T  ah* <*  v- id?  begin¬ 
ning  ."to  •  be  exposed  or-  ‘he  eroded,  surface  of . specimen  '>T,. 

WheP  th*»  material  overlying  ?  void  is  thin  enough,  .t*wr 
be  a  localized  weak  area  ur.;-h  will  erode- muon  *<-••>** 
than  material  volumes  with  no  por,***'.-.  However  the  total 
effect  for  the  number  of  l:-  .*  -*;.s  rr.t-  ■  M 

a  relatively  small  ' at.icr  r  •.  -  ,">rf  i.  »=rc-triu>  ra*-*-. 

The  approx1*  *  • *•  vrr=<-s<  -*?  !  ;e  •  -•.>**!.  ing^sfierial  for  the 
cr...;  :  i c’.cu  •  ir:  (U  is  mils. 
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Figure  12  Condition  of  eroded  surface  of 
at  the  end  of  test.  (1.75X) 


specimen  C-4 
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The  region  just  below  the  eroded  surface 
was  Investigated  in  greater  detail  in  order  to  asses  the 
extent  of  the  interaction  between  the  impulsive  loadings  on 
the  surface. and  voids  near  the  surface.  No  interaction 
effects  were  observed  in  the  cross  sections  examined.  No 
cracks  were  found  connecting  the  voids  with  the  surface  for 
the  sand  particle  impacts. 


5.0  GENERAL  DISCUSSION 

The  analysis  of  multiple  direct  impacts  of  glass  beads 
on  glass  plates  can  now  be  advanced  on  the  basis  of  the  micro¬ 
scopic  observations  and  the  general  form  of  the  weight-loss 
%  data.  It  was.  found  that  the  volume  of  material  removed  from 
the  specimen  at  any  time  t  .could  be  divided  into  three  com- 
.  ponents  as  indicated  in  Eq. (2)  v 

AV(t)  =  AV1(t)+AV2(t)+AV3(t)  . 

The  explicit  form  of  AV^(t)  was  derived  in  Section  3.2.2 

t 

AVx(t)  =  vdJ  K1(t-r)g(D)n<?(r)(l-AA1(T)-N(T)Ao)dT  (17) 

where  AV1(t)  is  the  volume  of  material  removed  per  unit  sur¬ 
face  area  of  the  specimen.  The  growth- rate  function  K^(t) 
denotes  the  rate  at  which  area  is  removed  once  a  pit  is  nucle¬ 
ated  and  the  sum  of  the  remaining  terms  in  Eq.(17)  represents 
the  number  of  pits  nucleated  per  unit  area.  Means  for  obtain¬ 
ing  the  parameters  appearing  in  Eq.(17)  were  described  in 
Section  3.2.3. 

As  the  exposure  of  the  specimen  continues  a  time  will 

eventually  be  reached  when  the  volume  of  the  specimen  assigned 

to  V, (t )  will  be  exhausted.  If  this  value  of  the  exposure 

* 

time  is  denoted  by  t^  ,  then 

t-t^V1)  3  riJj  .  (41) 

in  accordance  with  Eq.(12).  This  condition  requires  that  the 
Integral  in  Eq.(17)  equals  unity  when  t.«t^  .  We  further  note 
that  at  this  time 

AA1(t3<  )  ».l-N(t1t)A0  (4o) 
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where  N(t^*)  is  the  number  of  frustrums  protruding  from  the 
surface  of  the  specimen  per  unit  area  and  A^  is  the -aver age  v 
area  of  the  ring  fractures  on  the  specimen's  surface.  ? 

On  the  basis  of  the  experimental  observations,  volume  is 
not  removed  from  the  conical  frustrums  until  the  surrounding 
material  has  been  removed.  A  reasonable  ;f  orm  f Or  AYgt t )  is. 


The  lower  limit  in  the  integral  can  be  set  equal' to  zero, 
however  in  some  cases  it  may  be  advantageous  to  have  a  non- 
zero  value.  The  function  K2(p;)  is  the  time  rate  of  volume 
removed  per  nucleated  fracture ; resulting  in  material  removal 
from  an  exposed  conical  frustrum.  The  sum  of  the  remaining 
terms  in  Eq. (^3)  represents  the  number  of  critical  fractures 
nucleated  per  unit  area.  The  probability  function  <P2(t)  is 
time-dependent,  since  it  is  found  that  a  protruding  frustrum 
becomes  more  susceptible  to  decomposition  after  a  distribu¬ 
tion  of  fractures  develops  within  it  due  to  repeated  particle 
impacts.  A  second  time-dependent  factor  has  been  incorporated 
into  ?2(t )  which  is  the  fact  that  only  a  portion  of  the  area 
occupied  by  the  frustrums  receives  direct  impacts  from  the 
glass  bead3  due  to  the  geometry  of  the  frustrum  distribution 
as  shown  in  Fig.  6(c),  As  the  cones  disintegrate  an  increas¬ 
ing  area  is  exposed  to  direct  impacts.  An  appropriate  term 
for  the  function  ^2 (t)  would  be  a  sigmoidal  shaped  curve 
which  asymptotically  approaches  a  constant  value  >.*  which  Is 
always  less  than  or  equal  to  unity.  The  constant  value  of 
<Pg  is  reached  when  the  surface  of  the  specimen  is  highly  frac¬ 
tured  and  small  fragments  are  easily  removed.  At  an  exposure 
time  U*,  AAj  reaches  the  constant  value  indicated  in  Eq* (42) . 
Similarly,  there  will  be  a  time  when  <p0  is  also  constant,  then 
Eq.  ('43)  becomes 


AVp(t-.)  -  /  Kj?(D)n^(l-N(t1!‘)A0)dr 


=  Kg(D)n  vty*  (l-I^t^*  )  AQ  )  (t- t  ^  )=(oonst .  )  (t.-^  )  (44) 


or  the  volume  removed  asymptotically  approaches  a  linear  func¬ 
tion  of  time.  A  constant  rate  of  material  removal  is  found 
to  prevail  as  evidenced  in  the  form  of  the  weight-loss  data 
once  most  of  the  material  has  been  removed  from  the  surface 
layer  of  the  specimen.  As  the  exposure  to  the  erosive  envi¬ 
ronment.  continues,  a  point  is  reached  where  the  entire  surface 
layer  to  a  depth  T)  is  removed,  then 


t^AVt)  =  (1‘r>)D  (;+5) 

|  2 

In  the  derivation  of  Eq . (43;  we  have  assumed  that  all 
the  fragments  removed  from  the  surface  of  the  specimen  have 
the  same  volume  at  all  stages  of  the  frustrum  decomposition 
and  that  the  surface  area  of  the  specimen  which  participates 
in  the  erosion  process  can  be  represented  by  aA-j  ( t ) .  This 
latter  assumption  is  only  an  approximation  to  the  actual 
conditions,  however  it  becomes  quite  involved  to  include  an 
estimate  of  the  frustrums  included  in  the  area  AA-^  for  any 
generic  time  t.  This  leads  to  the  limiting  form  of  AV^(t) 

Ir.  Eq.(44)  where,  in  a  more  realistic  representation  (1-N(t^*) 
Ac)  should  be  replaced  by  unity. 

The  microscopic  investigations  do  not  indicate  that 
there  is  a  distinction  between  the  removal  of  material  from 
the  bulk  specimen  and  that  near  the  base  of  the  frustrums. 

This  observation  implies  that 

AV?  ( t )  =  J  K2(D)rn»2*AA2(r  )dr  (46) 

f2 

where  the  lower  limit  is  again  taken  to  be  a  non- zero  value 
signifying  the  exposure  time  when  this  final  phase  of  the 


-170- 


erosion  process  becomes  operative.  The  exposed  surface  area, 
AApft).  at  'he.  has-  t-r  the  surface  layer  is  difficult  to  des¬ 
cribe.  The  followin'1'  approximation  is  introduced  into  Eq . ( 46  ) , 


AA.,(t) 


av2(t) 

V^Tfn 


(^7 ) 


'■’hen  A Vp  =  AV^ft^'  ) ,  AAp(t)  =  1.  Hence  after  the  surface 
layer  has  been  removed 


AV3(t)  =  K2(D)n?2*(t..t2) 


(48) 


which  is  simply  a  linear  function  of  time  and  a  continuation 
of  the  latter  starr.es  of  the  erosion  process  contributing  to 
AV0 ( +  ) .  The  evaluation  of  all  the  parameters  appearing  in 
AV,  in  Eq. (46)  will  be  made  in  conjunction  with  Eq.(43). 


The  volume  removed  from  the  specimen,  or  equivalently 
the  weight-loss,  car.  be  computed  once  the  parameters  appear- 
lng  in  Eq.(17),  (43),  and  (46)  are  known.  The  required  data 
and  its  acquisition  for  evaluating  Eq.(17)  was  described  in 
Section  3.2.3.  The  additional  parameters  introduced  in  Eq.(43) 
and  (46)  are  Kp(D)  and<pp(t). 


The  expression  for  the  fracture  cone  penetration  depth 
D  riven  in  Eq.(?4)  shows  that  1  +  is  related  to  the  impact 
velocity,  the  radius  and  density  of  the  impacting  beads,  and 
the  elastic  constants  for  the  bead  and  target  materials 
through  the  Hertzian  theory  of  impact.  The  magnitude  of  P 
is  a  function  of  the  properties  of  the  specimen  in  its  initial 
state.  However  the  functions  Kp  and  <Pp  are  related  1:o  the 
highly- fractured  surface  of  the  specimen  during  the  final 
stage  of  the  erosion  process  as  shown  in  Fig.6(d).  It  is 
very  difficult  to  ascribe  mechanical  proper!. ies  to  the  highly- 
eroded  surface,  so  that  at  the  present  time  Kp  and  <P p  will 
have  !o  be  evaluated  experiment  ally .  Future  work  will  explore 
suitable  analytical  representations  for  these  functions. 
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Pho  above  model  1  r.  applicable  to  oblique  impacts  by 
glass  beads,  since  no  significant  differences  were  observed 
in  the  general  form  of  the  erosion  characteristics  for  this 
case.  The  weight-loss  data  indicated  that,  under  nearly 
equivalent  conditions  used  in  the  case  of  direct  impacts,  the 
magnitude  of  the  weight  loss  at  the  same  normal  velocity  and 
exposure  times  was  two  to  three  times  greater  for  the  oblique 
impact  as  compared  to  direct  impact.  The  ring  fracture  dia¬ 
meters  were  found  to  be  approximately  equal  for  both  impact 
conditions  and  in  general  agreement  with  the  values  of  the 
maximum  contact  diameters  calculated  from  the  Hertzian  theory 
of  impact,  however  the  penetration  depths  of  the  cone  frac¬ 
tures  were  noticeably  less  for  the  oblique  impacts  especially 
as  the  impact  velocity  increased.  It  is  concluded,  for  a 
reason  yr  t  to  be  discovered,  that  the  differences  in  the 
weight  loss  data  are  du '  to  a  greatly  increased  probability 
for  producing  ring  fractures  for  oblique  impacts.  However 
these  differences  can  be  readily  assimilated  in  the  general 
modeling  Context  developed  for  direct  impacts. 

The  direct  and  oblique  Impacts  of  irregularly- shaped 
solid  particles  dictates  that  more  extensive  modifications  of 
the  basic  analytical  model  are  required.  At  moderate  impact 
velocities  each  particle  is  capable  of  removing  material  from 
the  surface  of  the  target.  The  weight-loss  data  indicate  that 
a  steady-state  erosion  rate  is  established  after  a  very  short 
incubation  period.  The  irregularities  in  the  weight-loss 
data  for  both  direct  and  oblique  impacts  during  the  incuba¬ 
tion  period  require  additional  research  for  clarif i.eat ion . 

For  both  impact  conditions  the  normal  component  of  the  par¬ 
ticle  impact  velocity  was  held  constant,  and  it  was  found 
that  the  steady-state  erosion  rates  were  nearly  indentical. 

The  slightly  greater  values  in  Table  IX  for  the  oblique 
Impacts  may  be  due  to  the  fact  that  the  specimen  receives 
approx Irnat ely  30  percent  more  particle  impacts  per  revolution 
in  this  case  as  indicated  by  the  estimates  provided  in 


Table  V.  However  microscopic  examination  of  the  erosion 
damage  occurring  during  the  early  stages  of  the  erosion  pro¬ 
cess  did  not  reveal  any  significant  difference  between  the 
direct  and  oblique  impacts. 

Due  to  the  short  duration  of  the  incubation  period  and 
the  rapid  approach  to  a  uniform  rate  of  material  removal,  a 
general  relation  of  the  form  given  in  Eq. (48 )  is  an  appro¬ 
priate  description  of  the  weight-loss  data.  Of  course,  the 
fonn  of  the  coefficients  on  the  right-hand  side  of  Eq. (48  ) 
will  have  to  be  modified.  Some  guidance  in  this  regard  will 
be  obtained  from  the  explicit  form  developed  in  the  case  of 
the  70-micron  glass  beads.  However  the  characterization  of 
irregularly-shaped  particles  in  a  manner  suitable  for  incor¬ 
poration  into  Eq.  (48)  has  to  be  considered  if  this  relation 
is  to  have  general  applicability.  While  it  is  possible  to 
obtain  statistical  data  pertaining  to  the  size  of  irregularly- 
shaped  particles,  it  is  a  very  difficult  task  to  obtain 
appropriate  representations  of  particle  shape  as  indicated  in 
Section  3*4.1.  The  lack  of  an  adequate  representation  for  the 
particle  shape  was  pointed  out  as  a  contribution  to  the  failure 
of  the  correlation  between  the  erosion  rates  and  a  simple 
model  developed  by  Head  and  Harr.  This  does  not  exclude 

consideration  of  the  material  properties  of  the  irregularly- 
shaped  particles  which  will  also  have  a  significant  effect 
on  the  erosion  rates.  These  problems  require  additional 
attention. 

The  derivations  of  a  general  theory  for  the  erosion  of 
glass  plates  by  glass  beads  indicates  the  way  in  which  the 
material  properties,  particle  characteristics,  and  impact 
parameters  influence  the  material  removal  process.  While  glass 
beads  are  not  an  erosive  medium  commonly  encountered  in 
practice,  the  mathematical  analysis  of  the  impact  of  glass 
beads  on  the  pristine  surface  of  the  glass  plates  can  be 
given  in  considerable  detail.  The  more  highly  eroded  the 
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surface  of  a  brittle  material  becomes,  the  original  surface 
conditions  no  longer  prevail  and  the  problem  is  less  suscep¬ 
tible  to  mathematical  analysis.  However  the  value  of  the 
above  approach  is  that  it  provides  insight  into  the  general 
structure  for  a  predictive  erosion  theory  and  allows  every 
feature  of  the  basic  model  to  be  more  explicit  than  would  be 
possible  in  a  more  complex  erosion  process.  The  complexity 
of  the  material  response  and  erosive  medium  will  be  increased 
as  in  our  exploratory  investigations  of  sand  and  liquid  par¬ 
ticle  impacts  on  glass  plates. 


A  means  for  evaluating  the  pressure  distribution  imparted 
to  a  deformable  surface  by  a  spherical  liquid  drop  was  out¬ 
lined  in  Section  3.5.1.  The  transient  stress  distributions 
with  an  elastic  half-space  corresponding  to  the  applied  pres¬ 
sure  pulse  can  be  determined  in  an  approximate  way  using 
Blowers'  analysis.^2)  it  is  estimated,  on  the  basis  of  the 
limited  numerical  data  available,  that  the  intensity  ol  the 
disturbances  propagating  into  the  elastic  medium  due  to  a 
2mm  drop  impacting  the  elastic  body  at  velocities  up  to  1000  fp 
are  considerably  less  severe  than  the  results  obtained  from 
simple,  one-dimensional  analyses  would  indicate.  It  cannot 
be  established  that  the  transient  stress  distributions  due 
to  a  spherical  liquid  drop  impacting  a  deformable  surface  can 
be  represented  in  terms  of  the  equilibrium  stress  distributions 
shortly  after  the  drop  first  touches  the  surface  as  suggested 
by  Peterson^6)  until  after  additional  numerical  calculations 
employing  the  approach  described  in  Section  3*5.1  have  been 
completed  for  contact  times  in  the  vicinity  of  Igs.  The 
results  from  the  analysis  of  the  transient  stresses  due  to 
liquid  drop  impacts  on  deformable  materials  will  replace  the 
Hertzian  analysis,  if  necessary,  in  the  erosion  modeling 
studies,  so  the  rain  erosion  of  brittle  materials  can  be 
considered.  The  erosion  tests  established  that  material  remo¬ 
val  from  glass  plates  in  a  rain  environment  was  a  preferential 
process  of  pit  nueleation  and  growth,  therfore  the  erosion 
process  conforms  with  the  basic  premise  of  our  general  erosion 
model . 
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The  erosion  process  for  a  typical,  brittle  material  can 
be  considered  in  terms  of  two  major  subdivisions:  the  incu¬ 
bation  period  and  the  period  of  uniform  material  removal. 

The  erosion  resistance  of  a  specific  material  class  can  be 
improved  in  two  ways.  Once  the  material  parameters  affecting  . 
the  removal  of  material  in  the  incubation  period  are  identi¬ 
fied,  they  can  be  selected  so  the  duration  of  the  incubation 
period  is  extended.  This  means  that  the  magnitude  of  the 
weight  loss  at  a  particular  exposure  time  will  be  less  for 
the  modified  material  than  for  the  original  material.  The 
second  way  to  improve  the  erosion  resistance  is  to  select 
materials  which  have  material  properties  that  decreased  the 
magnitude  of  the  steady- state  erosion  rates.  If  the  opera¬ 
tive  erosion  mechanism  during  the  incubation  period  is  the 
same  as  that  in  the  period  when  the  steady-state  erosion  rate 
is  achieved,  then  modification  of  the  material  parameters 
associated  with  material  removal  in  the  incubation  period 
will  also  influence  the  steady-state  erosion  rates.  This 
would  be  true  for  the  solid  particle  erosion  of  brittle  and 
semi-brittle  materials.  Hence,  understanding  and  analyzing 
the  events  leading  up  to  the  initiation  of  material  removal 
Is  a  viable  approach  to  prescribing  the  material  character¬ 
istics  which  affect  erosion  behavior. 

The  erosion  mechanisms  have  been  identified  for  the 
TJHMW  polyethylene  and  the  graphite-fiber  reinforced  epoxy 
specimens  as  described  In  Section  4.1  and  4.3*  The  sand  ero¬ 
sion  of  UHMW  polyethylene  is  due  primarily  to  the  cutting 
action  of  the  irregularly- shaped  particles  striking  the  sur¬ 
face  of  the  aerodynamically- shaped  specimens  away  irom  the 
leading  edge.  The  direct  Impacts  of  the  sand  particles 
along  the  leading  edge  simply  produce  permanent  indentations 
in  the  polyethylene,  but  this  process  is  not  as  effective 
in  removing  material.  The  sand  erosion  of  the  graphite-fiber 
reinforced  epoxy  specimens  did  not  exhibit  a  significant 
angular  dependence;  material  removal  from  the  surface  was 


equally  effective  for  both  direct  impacts  and  cutting  defor¬ 
mations. 

The  rain  erosion  mechanisms  for  the  graphite-fiber  rein¬ 
forced  epoxy  are  particularly  interesting.  The  process  of 
material  removal  results  from  the  following  sequence  of  events. 

A  water  drop  striking  the  specimen's  surface  will  crack  the 
epoxy  surface  layer  through  direct  impact  or  upon  impacting 
in  one  of  the  pre-existing  shallow  depressions  on  the  surface. 
The  cracks  propagate  along  the  fiber/matrix  interface  when¬ 
ever  possible  which  is  an  area  of  weakness  in  the  composite 
structure.  Small  pitted  regions  appear  on  the  specimen's  sur¬ 
face.  Additional  cracks  exist  within  the  pits  and  are  enlarged 
through  the  pressures  applied  by  the  direct  impact  and  lateral 
outflow  from  drops  striking  in  the  vicinity  of  the  cracked 
area.  These  loadings  momentarily  force  the  two  faces  of  the 
crack  apart  and  it  extends  further  into  the  matrix  material. 

The  fracture  mode  in  both  the  initiation  stage  of  crack  develop¬ 
ment  and  the  penetration  of  cracks  into  the  matrix  corresponds 
to  mode  I  deformations.  Longitudinal  cracks  lying  parallel 
to  the  direction  of  the  carbon  fibers  are  readily  found  on 
the  eroded  surface.  Eventually  large  cracks  penetrate  the 
composite  isolating  small  bundles  of  fibers  and  matrix  which 
ultimately  are  only  attached  to  the  specimen  at  the  ends  of 
a  short  segment.  The  isolated  bundles  may  consist  of  a  few 
fibers  or  be  composed  of  a  grouping  which  is  five  to  ton  layers 
thick  and  an  equivalent  range  in  breadth.  Additional  droplet 
impacts  penetrate  the  severed  region  below  the  small  segment 
and  force  it  outward  from  the  surface.  At  this  point  the 
freed  segment  can  be  thought  of  as  a  rectangular  beam  fixed 
at  both  ends  and  loaded  by  a  distributed  load  on  one  of  its 
longitudinal  faces.  Removal  of  the  segment  results  through 
fracture  at  its  ends.  It  has  not  been  determined  whether  the 
bending  or  transverse  shear  stresses  represent  the  most  severe 
condi i ion. 
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The  microscopic  observation  of  1  ho  eroded  surface  am] 
the  above  reasoning  lead  to  two  important  aspects  of  the  mate¬ 
rial  removal  process.  ivLosi  *  fie  fracture  toughness  for  mode 
1  do  format  ions  of  the  epoxy  matrix  should  control  the  Initia¬ 
tion  of  the  fracturing  process.  Tt  is  only  the  mechanical 
behavior  of  the  matrix  material  and  not  the  mechanical  pro¬ 
perties  of  the  composite  which  govern  the  fracture  behavior 
under  these  conditions.  However,  it  must  be  remembered  that 
the  fracture  surface  tends  to  propagate  along  the  matrix/ 
fiber  interface,  so  the  bond  established  between  these  two 
components  of  the  composite  introduces  a  locally  weakened 
area  in  the  otherwise  homogeneous  matrix  material.  Secondly, 
the  initial  removal  of  small  segments  of  the  composite 
involves  ei’h.-r  mod,  .!  K,>  mode  TT  deformat. ions,  but  in  this 
case  the  crack  is  propagating  in  a  direction  transverse  1 o 
the  longitudinal  fiber  axes.  The  fracture  toughness  values 
must;  therefore  correspond  to  this  condition.  This  process 
continues  as  described  in  Section  4.3.1. 

Having  established  the  fundamental  mode  of  failure  of 
the  composite  at.  an  intermediat  e  to  advanced  stage  of  the 
erosion  process,  a  micromechanical  model  can  be  developed 
which  reflects  these  findings  in  quantitative  terms.  The 
first  approximation  provided  here  greatly  simplifies  the  act¬ 
ual  state  of  affairs.  A  number  of  conjectures  must  be  intro¬ 
duced  in  order  to  arrive  at  a  manageable  problem.  A  meaning¬ 
ful  analy* leal  approach  will  be  developed  in  future  research 
on  composites. 

The  nickel-coated  composite  specimens  exposed  to  a  sand 
environment  failed  by  a  loss  of  adhesion  at.  the  interface 
between  the  coating  and  Hie  substrate.  The  coal ing  then 
wrinkled  as  seen  in  Fig.  r>7  and  bH.  The  average  thickness  of 
the  coating  was  initially  3  mils  for  specimen  P-3-3  and  some¬ 
what  less  than  (>  mils  for  specimen  p-b-f>,  on  I  -p  basis  of 
f  he  weight  loss  for  '-.vih  specimen  at  Mi*-  time  the  coatings 
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buckled  and  i  he  area  of  the  specimen  exposed  to  solid  particle 
impacts,  the  uniform  reduction  in  the  coat  in/'  thickness  is 
only  oj  the  order  of  a  few  tenths  of  a  mil  which  would  hardly 
reduce  the  thickness  of  the  6  mil  coaling  to  a  critical  thick¬ 
ness  comparable  t:o  that  for  the  3  mil  coat  ini'.  The  buckled 
pattern  shown  in  Fig.  37  and  38  is  identical  to  the  buckling 
mode  for  a  thin  cylindrical  shell  of  ductile  metal  subjected 
to  uniaxial  compression.  This  observation  implies  that  there 
is  a  loss  of  shear  resistance  at  the  interface.  Tt  would  seem 
that  the  centrifugal  forces  acting  on  the  specimen  when  it  is 
whirling  at  the  end  of  the  rotating  arm  in  the  AFML-Pell  ero¬ 
sion  facility  contribute  to  the  failure  in  tin.  s  case.  All  hough 
our  calculations  showed  that  the  average  thickness  of  the  3 
and  8  mil  cnsi tug  did  not  reach  the  same  critical  value  before 
failure  occurred,  the  centrifugal  force  acting  on  the  o  mil 
coating  would  be  greater  than  nat  acting  on  1  he  3  mil  coat¬ 
ing.  The  critical  coating  thickness  would  therefore  have  to 
be  computed  accounting  for  this  difference.  Support  for  our 
conjecture  that  sand  erosion  reduces  the  nickel  coating  to 
a  critical  thickness  for  both  coating  thicknesses  is  obtained 
by  comparing,  the  different  failure  modes  found  in  the  rain 
erosion  tests.  Fig.  30  and  54,  whe/v  nc  measurable  weight  loss 
v/ as  recorded. 

The  research  conducted  within  this  program  led  to  several 
preliminary  concepts  for’  erosion  modeling  studies  in  a  variety 
of  materials.  The  proposals  made  here  will  be  developed  in 
future  research.  The  analytical  modeling  of  signif leant  fea¬ 
tures  of  the  operative  erosion  mechanisms  at  vari  ous  st  ages  of 
the  erosion  process  is  beginning  to  establish  the  connection 
between  material  properti  es  and  f  he  erosion  behavior  of  speci¬ 
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The  pxplir.Ii  Corn  for  a  complet  e  erosion  1  heory  has 
been  derived  for  cU red  and  oblique  impacts  of 
ill  ass  bead:',  >  r  glass  plates.  Future  research  will 
concent  raw  on  replacing  experimentally-determined 
parameters  appearing  in  this  theory  with  analytical 
expressions’-  containing  the  iinpac*  parameters  and 
material  properties  of  the  eroding,  particles  and 
target . 

Exploratory  investigations  were  initiated  for  the 
direct  and  oblique  impacts  of  sand  particles  on 
brittle  materials.  A  theory  of  erosion  will  he 
developed  for  these  conditions  based  on  the  experi¬ 
mental  observations  and  the  explicit  form  of  the 
theory  for  brittle  materials  eroded  by  small 
spherical  particles. 

Advancement  has  been  made  in  evaluating  the  transient 
stress  distribution  in  an  elastic  half-space  sub¬ 
jected  to  a  time-  and  spat ially-dopendent  pressure 
locally  applied  to  its  surface  and  approximating  a 
liquid  drop  impact.  The  results  obtained  from 
this  analysis  will  be  used  to  extend  the  erosion 
model,  to  liquid  drop  impacts. 

Frosion  tents  confirmed  that,  the  rain  erosion  of  glass 
plates  v/ as  a  pit  nucleat  ion  and  growth  process,  and 
therefore  glass  plates  provide  an  ideally  brittle 
response  for  extension  of  1  he  basic  erosion  model 
to  liquid  drop  impacts. 

T  e  erosion  mechanisms  for  HHMW  polyethylene  exposed 
to  sand  arid  rain  erosion  have  been  1  rienti.fi ed . 
Material  removal  by  sand  erosion  is  due  primarily 
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to  1  ho  out!  .in?  act  Ion  of*  the  par1  teles  striking 
the  targe?  obliquely.  Rain  erosion  of  UHMW  poly¬ 
ethylene  in  a  highly  preferential  process  of  pit 
nucleation  and  growth  which  conforms  to  the  basic 
hypothesis  used  in  the  crosi on  modeling  studies  of 
brittle  materials. 

The  erosion  of  nickel-coated  glass-fiber  reinforced 
epoxy  specimens  requires  further  investigation. 

It  is  suggested  that  an  interaction  occurs  between 
the  centrifugal  forces  applied  to  the  specimen  on 
the  rotating  arm  and  the  stability  of  the  reduced 
coating  thickness  as  material  is  eroded  from  the 
surface.  The  shear  resistance  at  the  coating  sub¬ 
strate  interface  appears  to  be  quite  low,  and  an 
analysis  based  on  the  buckling  of  a  thin  cylindri¬ 
cal  shell  in  compression  is  proposed.  A  different, 
form  of  failure,  cracking  and  localized  pitting  of 
the  coating,  was  found  for  rai  n  erosion.  The 
weight  loss  was  negligible  in  this  case,  and  the 
buckling  instability  did  not  appear.  In  general, 
the  stress  wave  interactions  at  the  coating/sub- 
strafe  interface  due  to  locally  applied  pressure 
pulses  on  the  surface  for  both  sand  particles  and 
rain  drops  should  be  considered  in  understanding 
adhesion  failures  for  coatings. 

The  erosion  mechanisms  for  graphit  '•-fiber  reinforced 
epoxy  specimens  exposed  to  sand  and  rain  erosion 
have  been  identified.  The  rat e  of  material  removal 
by  sand  is  the  same  due  to  direct  deformation  and 
the  cutting  action  of  i he  particle  striking  obliquely. 
In  +his  regard  the  nature  of  the  erosion  process 
is  the  same  as  that  for  glass  plates  eroded  by  sand 
part  icles.  The  analytical  description  of  the  In*  tor 
case  will  therefore  be  direoily  applicable  ) o  ■ he 


-180- 


sand  '-rosJon  of  the  compos]  *  '>  material  which  behaves 
in  a  i  r ■  •  i  1  o  manner.  Th--  rain  erosion  mechanisms 
were  quite  different  from  the  mechanisms  associated 
with  sand.  A  significant  material  property  which 
affects  the  •-rep ion  behavior  is  the  fracture  tough¬ 
ness  of  1  lie  epoxy  matrix  and  the  composite  material. 
A  ml  eromeebanical  model  is  proposed  for  relating, 
fracture  toughness  to  the  erosion  process. 
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